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Development of the Lockheed P-80A 
Jet Fighter Airplane 


CLARENCE L. JOHNSON* 
Lockheed Aircraft Corporation 


INTRODUCTION 


~~ PURPOSE OF THIS PAPER is to present a summary 
of the development of the Lockheed P-80 ‘‘Shoot- 
ing Star’ jet-propelled fighter airplane, with particular 
reference to the interesting engineering problems en- 
countered and the means used to obtain their solutions. 
As this airplane is the first American jet-propelled tac- 
tical aircraft, many lessons have been learned from its 
design, testing, and operation which should be of great 
assistance in improving future similar aircraft. Because 
of the great secrecy surrounding this project throughout 
the war, it was not possible to discuss many of the 
interesting features of the P-80 at an earlier date. 


Background for Type 


The Army Air Forces assigned Lockheed the task of 
building a jet-propelled fighter on June 23, 1943, when 
it became apparent that only an airplane powered by 
a turbojet power plant could compete with the new 
German types such as the ME 262. Lockheed had 
demonstrated interest in jet aircraft since 1939 by ac- 
tual design work on a radical fighter equipped with an 
engine of its own conception. Time did not exist, how- 
ever, to build this type, and the gamble did not seem 
warranted-to the Army Air Forces under wartime pres- 
sure. It was, therefore, decided to construct a more 
conventional type, utilizing first an engine of English 
manufacture. 


The Prototypes 


Under conditions of the utmost cooperation from the 
Army Air Forces, Lockheed designed, built, and com- 
pleted preliminary ground tests on the XP-80 aircraft 


Presented at a closed session, Fifteenth Annual Meeting, 
LA.S., New York, January 28-30, 1947. Received July 25, 1947. 
* Chief Research Engineer. 


in 141 days. An accident to the intake ducts, resulting 
in damage to the original engine, prevented flying 
for several more weeks because of the lack of a power 
plant. 

The XP-80 (Fig. 1) utilized successfully a British de 
Havilland Halford engine having a 2,500-lb. static 


Fic. 1. Prototype XP-80 with Halford jet engine. 


Fic. 2. Prototype XP-80A with General Electric I-40 engine. 
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Fic. 3. View of detachable fuselage section sliding over engine. 


thrust rating. Production of the Halford engine was 
not forthcoming in time, however, to allow immediate 
procurement of P-80 type aircraft based on the original 
design. 

' The lack of Halford engines was overcome by the 
production of the General Electric I-40 type in this 
country. To use this engine, however, it was found 
necessary to redesign completely the XP-80 to make a 
new model known as the XP-80A. This airplane had a 
new fuselage of greater length, 25 per cent greater de- 
sign weight, new ducts, engine installation, gear, etc. 

The XP-80A (Fig. 2) was then designed, built, and 
flown in 138 days. Production followed closely behind 
the experimental articles, but, while a few aircraft got 
to the European theater of operations during the war, 
none actually engaged in combat. 


New Design Features 


The P-80A type airplane incorporates many new de- 
sign precedents (Figs. 3 through 11) which will no doubt 


affect later types. Some of these are: (1) detachable 


fuselage tail assembly, (2) quick-change engine installa- 
tion, (3) control surface hydraulic booster, (4) com- 
bined cockpit pressurization and cooling, (5) wing-tip 
fuel tank and bomb carrying arrangement, (6) fuselage 
dive brakes, (7) automatic fuel transfer system, (8) 
short chord control surfaces, (9) precision contour con- 
trol and special finish, (10) high acceleration type arma- 
ment installation, (11) jet engine water-injection sys- 
tem, (12) internal radio antennas, and (13) automatic 
emergency fuel system for engine. 

Problems encountered in developing the above and 
other features are discussed below. 

Fig. 5 shows the inboard profile of the airplane. It is 
largely self-explanatory. 


POWER-PLANT INSTALLATION 


Accessibility 


As in most new aircraft in the high-performance cate- 
gory, the power-plant installation presented the greatest 
number of development problems. The General Elec- 
tric I-40 engine is a double entry, radial compressor jet 


engine requiring a plenum chamber type of installation. 
This makes for a clean installation mechanically, and 
immediate steps were taken to provide a quick-change 
engine installation by designing the aft fuselage and tail 
assembly to be removed as a unit. Fig. 3 shows the 
parting line where three bolts are used to hold the tail 
section in place. Control cables have quick disconnects, 
as does the engine tailpipe. The engine may be run 
with the aft fuselage removed for inspection and test. 

In field use, Army crews have demonstrated engine 
changes in less than 20 min. The value of such a design 
feature to the armed forces cannot be overestimated, 
particularly when the experimental nature of the engine 
is considered. 


Air Duct System 


The lack of duct connections to the compressors does 
have disadvantages in performance, however. A loss 
of inlet ram, preheating of the intake air, and similar 
problems result in loss of net thrust. The question has 
often been asked as to why a nose inlet was not used 
for the ducts on the P-80A. The answers are these: 

(1) An inlet ram average of almost 100 per cent 
(Fig. 14) is obtained at the face of the present ducts. 
This cannot be bettered by a nose inlet. 

(2) The use of ducts close to the wing improves the 
critical Mach Number of the wing-fuselage intersection. 

(3) The armament installation was considerably 
improved with no duct interference. 

(4) Side inlets saved considerable weight. 

(5) A better fuselage nose shape for visibility is ob- 
tained. 


OUTLET 
(LOW TEMP) 


Fic. 4. Cockpit cooling turbine unit and aftercooler developed in 
joint tests by AiResearch and Lockheed for P-80A. 
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. Adjustable Light 

. Oxygen Cylinder 

. Ammunition Box 

. Armament Junction Box 


Command Radio 


. Instrument Panel 
. Bullet-Proof Windshield Panel 
. Gun Sight 


Seat 


. Fuel Level Gage 


FUSELAGE CONTENTS 


11. Fuel Tank 

12. Intake-Air Duct 

13. Command Radio Antenna 
14. Fuselage Aft-Section 


Attaching Joint 


15. Elevator Control Differential 
16. Air Speed Pitot 

17. Tail Pipe Support Track 

18. Tail Pipe 

19. Rernote Compass Transmitter 


20. Tail Pipe Clamp 

21. Elevator Tab Motor 
22. Engine 

23. Intake-Air Seal 

24. Engine Mounts 

25. Aileron Booster Unit 
26. Wing Spars 

27. Aileron Torque Tube 
28. Elevator Push-Pull Tube 
29. Identification Radio 


Sub-Cockpit Junction Box 


. Battery 


32. Identification Radio Antenna 


. Elevator and Aileron Control 
Assembly 


34. Nose Alighting Gear 
35. Rudder Pedals 
36. Fuselage Nose-Section 


37 
38 


Attaching Joint 
. Cartridge Case Ejection Doors 
. .50 Calbr Machine Guns (6) 


tion. 


Fic. 5. Inboard profile showing fuselage contents. 


Fic. 6. Flight photograph showing wing-tip fuel-tank installa- 


A number of problems did develop with the side inlets 
used. These were: 
Unstable duct flow at low flow ratios due to flow 
breakdown from boundary-layer separation. 


This 


Fic. 7. 


View of extended fuselage dive brakes. 


bd 


gave conditions of duct “rumble” and directional 
“snaking’’ until cured by boundary-layer bleeds. 
(2) Structural interferences necessitated curved 
ducts having poor expansion angles, which resulted in a 
ram loss greater than desired. This problem would not 
be solved using nose inlets. 
(3) Dumping losses in the plenum chamber, and ; 
particularly losses flowing over the engine to get to the 


rear compressor inlet, resulted in ram losses. 
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Fic. 11. View of left side of cockpit. 


Fic. 8. General view of tail surfaces showing use of narrow 
chord surfaces. 


INTAKE SYSTEM 


Fic. 12. P-80A air duct system showing boundary-layer bleed 
arrangement. 


Fic. 9. View showing nose armament. Note ammunition boxes 
over guns to get gravity feed. 


WITH BOUNDARY LAYER BLEED 
AND HORIZONTAL DUCT VANE. 


PLENUM CHAMBER PRESSURE 
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FROM ; SCALE MODEL LOCKHEED TUNNEL 
Fic. 10. View showing simplicity of cockpit arrangement. Fic. 13. 
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It seems inevitable that the use of axial flow compres- 
sors will result in simpler duct problems and give higher 
ram efficiencies than those obtainable on double entry 
centrifugal jet engines. 

The net ram finaliy obtained on the P-80A at high 
speed is 70 per cent of the available adiabatic compres- 
sion pressure. Figs. 12, 13, and 14 show the duct de- 
sign, the variation of ram with flow, and the benefit 
derived from the boundary-layer bleed close to the 
fuselage. Without the latter unit, rather amazing 
noises, buffeting, and directional stability were ob- 
tained. When the throttle was closed at high speed, 
the airflow into the ducts was suddenly reduced and a 
thick boundary layer developed on the inside of the 
duct at the fuselage wall. Uneven flow developed in 
the two ducts, sometimes even reversing in one. An 
organ-pipe tone resulted from the fluctuating flow which 
could be heard clearly on the ground when the airplane 
was a mile high. Stalling of the duct lips occurred 
which affected the flow over the fillets. A directional 
oscillation of 1'/2° and of high frequency (about 2 
sec.) known as “‘snaking’’ resulted. Strangely enough, 
“snaking’’ has developed on several jet aircraft, appar- 
ently from causes other than duct troubles. Generally, 
the fuselage shapes used on these aircraft contribute 
to the problem. 

Pressure surveys in flight and on a one-third scale 
wind-tunnel model resulted in the development of a 
boundary-layer bleed which by-passed the low energy 
air above and below the wing. Fig. 14 shows the im- 
provement in flow conditions obtained at the duct en- 
trances, and Fig. 13 shows the improved ram obtained 

especially at low flow ratios. This is important in that 
the air-starting characteristics were greatly improved 
as the windmilling engine r.p.m. doubled with the 
higher ram. Both duct rumble and ‘‘snaking’’ were 
eliminated with the installation of the boundary-layer 
bleed. A reduction in speed of 2 to 3 m.p.h. was a cheap 
price to pay for these gains. 

With a plenum chamber installation, consideration 
must be given to the suction forces developed on the 
ground with the engine running, since these are high 
enough to collapse ducts (as on the XP-80) or the fuse- 
lage. To prevent this, as well as to lower the negative 
ram, spring loaded doors (Fig. 15) were installed which 
prevented any appreciable suction from developing in 
the plenum chamber. They were entirely successful 
and do not flutter. 

As a further duct development, a flush-type installa- 
tion (Figs. 16 and 17) is being developed which should 
have an extremely high critical Mach Number as well as 
low drag. 


Tailpipe Problems 


The tailpipe installation presented numerous de- 
velopment problems. Considering the great tailpipe 
length (71/2 ft.), the low loss of 1 per cent thrust is 
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Plenum chamber pressure relief doors in fuselage be- 
hind canopy. 


Fic. 16. Front view of flush air ducts. 


Fic. 17. Side view of flush air ducts. 
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considered satisfactory. Tailpipe insulation inside the 
dura] fuselage required care, but several solutions de- 
veloped. The variation of production engines in thrust 
(up to 500 Ibs.) gave a most difficult problem, in that it 
became impossible to choose a fixed tailpipe nozzle 
that would fit all engines. Fig. 18 shows the extreme 
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importance of tailpipe nozzle size on engine thrust and 
tailpipe temperature. A decrease in nozzle diameter 
from 19 to 18'/2 in. increased the static thrust 400 Ibs. 
and the jet temperature 150°F. The change of airplane 
performance was naturally marked. A variable tailpipe 
nozzle area would, of course, be an easy solution to this 
problem, but unfortunately it is not available on the 
I-40 engine. Production tolerances on jet engines have 
a much greater effect on performance than they do on 
piston-type reciprocating engines. 

To investigate an early accident, a test was made ona 
complete airplane to determine whether fuselage deflec- 
tion in maneuvers could load up the tailpipe or engine 
to cause failure. Fig. 21 shows the aircraft on a jig 
which imposed a horizontal tail load of 8,000 Ibs., 
severe vibration to simulate buffeting, and which allowed 
the engine to be run at take-off thrust. In spite of 2 in, 
of fuselage deflection under this load, no damage was 
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Airplane set up for tailpipe load test with engine 
running. 
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obtained. A test was also run to see what tension load 
could be sustained by the tailpipe disconnect joint. 
Computed loads of several thousand pounds were only 
30 per cent of the joint capacity, but in spite of this a 
joint was developed which could resist 13,600 Ibs. of 
tension. 

To investigate what would happen with a failed baf- 
fle seal on the engine, particularly if the pilot closed his 
throttle at high speed, another series of tests was run. 
In the case noted, high ram pressures would be present 
outside the tailpipe and static pressure inside. A thin- 
gage pipe of large diameter can take only a small com- 
pression load. Fig. 22 shows such a failure at 7 in. 
mercury pressure differential. The tailpipe was con- 
structed to take the loads developed by this condi- 
tion. 


Engine Fuel System 


Present jet engines are equipped with a single fuel 
system that is perhaps the most vulnerable part of the 
whole engine with regard to reliability. For years, re- 
ciprocating engines have had dual ignition systems for 
added safety. Early in the development program of the 
P-80A airplane it became apparent that an emergency 
fuel system should be developed for the I-40 engine 
which could be used particularly on take-off and landing 
in event of failure of the larger elements of the standard 
engine system. This emergency fuel system is shown 
schematically on Fig. 23. An electrically driven high- 
pressure fuel pump is hooked into the engine fuel system 
in such a way that its fuel supply by-passes the three 
critical engine units which are subject to failure. A 
pressure valve is inserted in the system so that any 
time the fuel pressure of the main engine system falls 
below a predetermined value, the emergency fuel system 
provides sufficient fuel to maintain at least 90 per cent 


Fic. 22. Tailpipe failure under negative pressure differential 
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of take-off thrust. This emergency fuel system has 
saved a large number of airplanes, and as a result of its 
operation it is apparent that all jet and turbine type 
engines should be equipped with completely inde- 
pendent dual fuel systems to provide equivalent safety 
to that which now exists on reciprocating engines that 
have dual ignition. 

The writer has gone into considerable detail on the 
foregoing problems with the view that it might call the 
attention of others to the safety aspects involved which 
can readily be overlooked in new designs. 

During the testing period of the airplane it became 
increasingly apparent that the pilot should have some 
indication of the thrust developed by the engine during 
a run-up prior to take-off. With conventional engines 
and propellers, a check to take-off r.p.m. and manifold 
pressure in low pitch gives assurance that take-off 
power is available. With a jet engine there was no 
easy means for checking the same characteristic. A jet 
engine loses thrust approximately three times as fast 
as a reciprocating engine with increase in temperature. 
There was also a considerable loss in thrust with alti- 
tude. In an effort to find an easy means to instruct the 
pilot regarding available thrust, it was noticed that the 
burner ring fuel pressure was an excellent indication 
of thrust regardless of outside air temperature, normal 
altitude variations, and actual tailpipe nozzle size. 
Fig. 19 shows data obtained on a given engine when 
tested under varying conditions. If it were not for the 
basic variation of production engines, it would be easy 
to set minimum burner ring pressures for take-off 
which would assure the pilot ample thrust under all 
conditions. This may still be done if a different value 
is chosen for each particular engine. 


Airplane Fuel System 


The basic airplane fuel system consists of a number 
of wing tanks having an automatically operated transfer 
system to pump fuel to the main fuselage tank. A 
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manual override and a light system to indicate empty 
tanks are also available for emergency operation. 

The wing-tip fuel tanks are emptied by compressor 
air pressure. 


Fluid Injection Development 


In order to improve the take-off and military perform- 
ance of the P-80A airplane, Lockheed undertook to de- 
velop a water injection system for the I-40 engine. 
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The effect of water injection had been studied by other 
groups on other engines, but it had not been applied 
practically to any complete installation to the knowl- 
edge of the writer. It was, therefore, necessary to 
undertake a considerable research program on the ef- 
fect of water and water-alcohol injection into the engine 
compressor inlets, as well as to develop the mechanical 
means for assuring good distribution and proper flow 
rates. The effect of injection of a water spray into a 
jet engine is made up of several factors. Fig. 25 shows 
some of these in detail. Because of evaporation effects 
at the compressor inlet, the inlet air temperature is de- 
creased between 20° and 30°F. for high rates of injec- 
tion. This results in an increase in mass flow measured 
to be approximately 10 per cent by a tailpipe rake sur- 
vey. The evaporation of water during the compression 
cycle reduces the compressor outlet temperature con- 
siderably and necessitates the injection of substantially 
higher amounts of fuel to the combustion chambers, 
The mass of water injected turns to superheated steam 
in the tailpipe and tends to increase outlet velocity, 
Substantial thrust increases were developed, approach- 
ing 30 per cent of the static thrust for flow rates around 
25 gal. per min. for water. This flow rate was between 
two and two and one-half times the rate of fuel injec- 
tion. There has been some indication that the per cent 
increase in thrust is not a constant with varying ambi- 
ent air temperatures, in that higher rates of increase are 
obtained at high outside air temperatures and lower 
values are obtained with cold temperatures. The use 
of water injection reduced the distance required to 
clear a 50-ft. obstacle to 70 per cent of the distance re- 
quired with no fluid injection. 

A number of problems developed immediately when 
an attempt was made to use water injection on take-off 
and at altitude. On the P-80A, compressor pressure is 
used to drive the flight instruments, to pump fuel from 
the wing-tip tanks to the main fuel tank, and to pres- 
surize the cabin. Naturally, with water injection tak- 
ing place, the instruments become loaded with water, a 
considerable amount gets into the wing-tip fuel tanks, 
and the cockpit is filled with water. vapors causing 
windshield fogging. Furthermore, when water is 
carried in the airplane at altitudes above 15,000 ft., it 
will freeze solid because of the low temperatures. To 
overcome some of these difficulties a combination of 
water and alcohol was tested. The results obtained in- 
dicate that the use of alcohol reduces the effectiveness 
of the mixture in increasing the total thrust. Fig. 26 
shows this effect. While the presence of 50 per cent 
alcohol overcame the freezing problem, there still exists 
the problem of water getting into the wing-tip tanks 
and instruments, and the pilot problem is increased be- 
cause he now becomes drunk as well as blind. It is, 
therefore, necessary to shut off all air bleeds from the 
compressor when the water-alcohol combination is 
being used. Additional work is being carried on to ob- 
tain a solution to this problem. 
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Fuel Type 


The problem of using gasoline instead of kerosene 
for fuel constantly cropped up during the aircraft de- 
velopment. The advantages of kerosene are decreased 
fire hazard on the ground (this is questionable at alti- 
tude) and improvement of airplane range because kero- 
sene has more B.t.u.’s per cu.ft. than gasoline, so more 
fuel can be placed in a given space. Present jet engines 
with their high-pressure fuel systems and leaky joints 
in the plumbing, tailcone, and tailpipe offer consider- 
able hazard, particularly in ground starting, stopping, 
and taxiing. Kerosene is comparatively innocuous 
under this condition compared to gasoline. From an 
overall safety point of view the advantages then go to 
kerosene, although it is probable that gasoline would 
give improved starting and high-altitude burning char- 
acteristics. 


Starting 


The power required to start jet engines is consider- 
ably more than that required for reciprocating engines. 
Current surges on starting up to 900 amp. at 24 volts 
must be supplied for a short period during the starting 
cycle on the I-40 engine. This requires the use of a 
special starting cart, although batteries have been pro- 
vided by the Air Forces which are capable of giving 
about six to nine starts. Starting in the air must be 
done with caution, since it is easy to flood the tailpipe 
with fuel when making a false start at altitude. This 
presents a considerable fire hazard. Proper procedures 
for air starts and studies of the combustion problem 
at altitude required a considerable amount of flight 
time during the development period. 

The problem of turbine wheel failures and engine 
overspeeding will be referred to in another part of this 


paper. 
AERODYNAMIC DEVELOPMENT 


The aerodynamic development of the P-80A entailed 
the use of all types of wind tunnels and a considerable 
amount of flight research. Only the more interesting 
problems encountered will be discussed. 


Compressibility Effects 


A great deal of compressibility research has been 
done on the P-80A by Lockheed, and both the N.A.C.A. 
and the Army Air Forces are continuing to use the type 
for basic research. The effect of Mach Number on the 
profile drag variation of the P-80A is shown in Fig. 27. 
These data were derived in flight using tailpipe thrust 
and dive angle measurements, as well as horizontal ac- 
celerometers in the manner developed by the N.A.C.A. 
for deriving drag coefficients in dives. It is interesting 
to note that an increase of thrust of over 100 per cent is 
required to obtain a Mach Number of 0.80 compared to 
0.70. 
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Fig. 28 shows the data presented in a different form 
to illustrate the effect of ambient air temperature on 
level flight speed. 

The assumption is becoming current that jet airplanes 
can develop higher maximum level speeds with high air 
temperatures than at lower temperatures. This is due 
to the increase in true speed for a given Mach Number, 
assuming the aircraft reaches a constant Mach Number. 
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Point A on Fig. 28 indicates the high speed (605 
m.p.h.) obtainable with a given aircraft and engine on 
a standard day. The same combination on a 100°F.- 
temperature day loses in thrust but still can reach ap- 
proximately the same Mach Number, so it obtains a 
speed of 610 m.p.h. in spite of the power loss. 

Assuming the same airplane but a lower thrust engine, 
a high speed of 573 m.p.h. (Point C) can be reached on 
a 60°F. day. As the temperature increases to 100°F., 
speed is lost (567 m.p.h. at Point D) because of the 
variation in shape of the drag curve as the compressi- 
bility effect reduces. The power loss is predominant 
and a lower Mach Number is reached because of this 
limitation. 

It is therefore possible to have a gain, loss, or no 
change in true speed as the air temperature varies, de- 
pending on the shape of the drag and thrust curves for 
the particular aircraft-engine combination involved. 

During dive tests on the P-80A, many interesting 
phenomena were noted. Early in the test program it 
became obvious that manufacturing differences be- 
tween airplanes; wing roughness, particularly at the 
leading edge; and the rate at which speed was in- 
creased to the critical Mach Number, all had important 
bearing on the airplane response. Some airplanes 
tested developed a decided rolling tendency before 
showing any longitudinal pitching effect. This was 
probably due to unsymmetrical formation of wing 
shock waves. Others would go to higher Mach Num- 
bers and develop a gradual nose-down pitching moment. 
During recovery, high accelerations could be easily ob- 
tained when this diving moment was lost. One of the 
most remarkable effects noted was the great decrease 
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in critical Mach Number on an airplane that had en- 
countered heavy rain at high speed which severely pit- 
ted the wing leading-edge paint (see Fig. 55). This 
airplane encountered a strong diving tendency at a 
Mach Number 0.10 lower than that common for the 
type. 

During a flight with wing-tip fuel tanks on, the 
writer could see plainly the formation of shock waves 
at the wing-tank intersection under the existing at- 
mospheric conditions. Major Barsodi of the A.A.F. 
had previously noted similar shock waves on a P-51 in 
adive. Further study of the conditions under which the 
actual shock waves may be seen in flight might lead 
to a useful method for studying compressibility 
phenomena. 

Contrary to statements made in a leading national 
magazine, no compressibility effects have ever been 
noted when firing the nose guns on the P-80A, even at 
speeds over 600 m.p.h., nor have any mysterious vibra- 
tions been set up which caused the total disintegration 
of the aircraft. 


Wing-Tip Fuel Tanks 


The unusual location of the droppable wing-tip fuel 
tanks was the result of complete wind-tunnel tests on 
eight types of installations, having equal volumes, to 
determine the optimum aerodynamic and structural 
position. Problems of ground clearance and flap or 
aileron clearance and the requirement for carrying large 
1,000-lb. bombs as alternative load complicated the 
solution. The final underhanging wing-tip position 
was chosen for the following reasons: 

(1) The drag increase at high speed with tip tanks 
on was only 12 per cent of the airplane profile drag 
compared to 27 per cent for the P-38 pylon-type instal- 
lation or 30 per cent for the slipper type (see Fig. 29). 

(2) The wing-tip position increased the effective 
aspect ratio from 5.50 to 9.9 by providing a useful end 
plate on the wing. Note Fig. 29 where above a C, = 
0.45 the profile drag is Jess with the tanks on than with 
them off when based on the same aspect ratio. Actual 
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glide tests on the airplane substantiate this condition 
in that the airplane glides 92 miles power-off from 40,000 
ft. with tanks on and 87 miles in the same conditions, 
tanks off. 

(3) The tip location had minimum adverse effect on 
the airplane critical Mach Number. The P-80A has 
flown to speeds over M=0.75 with the tiptankson. In- 
board locations adversely affect the wing critical speed. 

(4) No dropping clearance problem (from tail or 
the wing itself) existed at the wing tip. 

(5) With the exception of local wing-tip regions, the 
structural strength of the airplane is increased for all 
flight conditions over the condition with tanks empty 
and off. This is due to the favorable dead weight ef- 
fect of the fuel. 

(6) The aileron effectiveness is increased consider- 
ably by the presence of the tip tanks. This compen- 
sates fully for the increased rotational moment of in- 
ertia. 

(7) If extreme rates of roll would inadvertently be 
developed, the tip tanks would throw clear with no 
damage to the aircraft. 

(8) While a centered position on the wing tip gave 
better dihedral effect than the low position (Fig. 30) 
and equal drag values, the inability to carry existing 
bombs in this manner and structural problems with the 
tanks themselves outweighed this one advantage. 

Perhaps the greatest disadvantage of the under- 
hanging tip tank position is the adverse effect on the ef- 
fective wing dihedral. Fig. 30 shows the results, which 
are confirmed by flight tests. While the centered tank 
location increases the dihedral, the underhanging posi- 
tion gives neutral stability. This condition was ac- 
cepted, however, to gain the advantages noted in Item 
8 above. 

During bombing tests, the Air Forces encountered a 
totally unexpected result from wing-tip vortexes. Cer- 
tain bombs had their tail fins unscrewed by the strong 
tip vortexes. When the tails were locked in place by 
pins through the threads, they were then broken off by 
the high loads encountered. 

The question is often asked about landing with a full 
tip fuel tank on one wing. Because of the powerful 
lateral control available through the hydraulic aileron 
booster, such a landing can be made safely. In the 
course of testing at Lockheed two such landings were 
made. It is not good judgment to do so, however, be- 
cause of the difficulty of holding the aircraft straight 
after landing and the possibility of dragging the tank 
and setting the fuel on fire. Many landings have been 
made with both wing-tip tanks on and full of fuel. 

For a long-distance record flight, one P-80A was fitted 
with two 330-gal. tip tanks. In spite of carrying 2,500 
lbs. on each tip, the airplane flew well. 


Aerodynamic Smoothness 


An extreme effort was made during the P-80A de- 
velopment to obtain aerodynamic smoothness over the 
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whole airplane. When it is considered that the total 
drag of the subject aircraft is only between 60 and 70 
per cent of that of a comparable internal-combustion- 
powered, propeller-driven fighter, the importance of 
smoothness is apparent. With no propeller to develop 
a turbulent airflow over the wing, every opportunity 
exists to take full advantage of laminar-flow airfoil de- 
sign. 

It was not expected that much gain would be ob- 
tained in high speed because of extreme smoothness, 
since it was apparent that the formation of shock waves 
would cause transition in any case. As long as reason- 
able leading-edge smoothness existed so thatlocalshocks 
were not precipitated, a high degree of smoothness could 
not be justified from the point of view of maximum 
speed. This was not true for speeds under the critical, 
however. A noticeable benefit could be obtained in 
range, cruising speed, and climb. 

To obtain these gains, Lockheed set up rigid manu- 
facturing tolerances and used the best available putties, 
fillers, and paints to get smoothness. A research pro- 
gram was carried on for two years to improve the qual- 
ity of finishes. A careful series of flight tests was made 
to evaluate surface finish. An airplane built carefully 
to eliminate poor skin joints and contours was flown with 
a selected jet engine. Speed and climb data were de- 
rived. The engine thrust was measured before, during, 
and after all tests. The airplane was then given a good 
finish utilizing putty, filler, and paint, and the tests 
were repeated. The data on speeds are shown in Fig. 
31. Correction was made to account for the gradual 
thrust decrease of 14 lbs. per hour. The results indicate 
about 7 m.p.h. gain in high speed at all altitudes and 
bear out the expectation that, when compressibility ef- 
fects on the wing limit speed, little is gained by extreme 

‘smoothness. At cruising r.p.m. (10,500 and less) the 
gain is marked, as a speed increase of 10 to 25 m.p.h. 
is obtained. 

These tests were run on a carefully constructed air- 
plane with a good paint job. Naturally, different results 
would be obtained with rougher aircraft, but the data 
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EAFECT PAINT ONT SPEED 
TION WITH ALTITUDE | 
AIRPLA 


4) “ENGINE 


FINISH + M. 


indicated are believed to be a good average for the type 
considered. Unfortunately, if the paint deteriorates 
unduly in service, the gains indicated are easily lost. 
It is probably in order to point out some of the require- 
ments for high-speed aircraft paints and fillers resulting 
from P-80A tests. These are: 

(1) Ability to withstand temperature changes from 
+160° to —80°F. in 3 min. or less. These tempera- 
tures are encountered in a zoom climb from sea level to 
35,000 ft. or, more rapidly, in a dive descent from 40,000 
to sea level. 

(2) Ability to maintain a firm bond to dural under 
above temperatures when subject to vibration and a 
large range of stress changes in the dural. 

(3) Ability to withstand severe rain at speeds of at 
least to 600 m.p.h. 

(4) Reasonable application means and time. 


RATE oF 
INDICATED, SPBED 


IN INS: 


Fic. 31. 


4 


RATE OF ROLL +DEGREES/SE: 
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In the course of testing available paints, some inter- 
esting data were obtained. Using a nozzle arrangement 
spraying oil to simulate encountering rain at 700 m.p.h., 
it was found that Alclad would erode badly and dural 
also would become rough. Metal spray directly over 
joints was not satisfactory as a practical means for ob- 
taining smoothness. 

Until better materials are available from the paint 
manufacturers, the results of Lockheed investigations 
indicate that true aerodynamic smoothness can be ob- 
tained only by using heavy plate construction and 
grinding to contour. 


Problems of Control 


No particularly difficult problems of control were en- 
countered within the normal Mach Number range of 
the P-80A. Extremely good maneuverability was ob- 
tained using a hydraulic booster of the type previously 
developed for the P-38 on the ailerons. Using a boost 
ratio of 15:1 (the mechanism multiplies the pilot effort 
at 15), a rate of roll over 150° per sec. can be reached at 
intermediate altitudes (Fig. 32) and much higher rates 
at 30,000 ft. Boost ratios of 28:1 have been flown suc- 
cessfully and phenomenal rates of roll obtained, but 
these were not considered required or desirable. The 
aileron cables are metal-clad equivalent to */s in. basic 
diameter. In spite of the precautions taken to reduce 
system deflection, only 50 per cent total aileron deflec- 
tion is obtained at high indicated speeds. 

No aerodynamic balance is used on any of the control 
surfaces. A variable mechanical advantage and spring 
tabs are used on the elevator to get desirable forces. 

No power effect due to the jet is apparent either longi- 
tudinally or directionally. 

Complete spin tests have been made on the P-80A. 
With the exception of the case with full wing-tip fuel 
tanks, recovery is excellent. In view of the fact that 
the tanks are easily jettisoned, no problem exists in 
this configuration. The spins developed are unusual 
in that the nose rises and falls considerably but the re- 
covery is easy and rapid, generally within one-half of a 
turn. No changes had to be made to the airplane to 

pass spin tests except to reinforce the elevator and rud- 
der counterweights to take over 150g. 


Fuselage Dive Flap 


To permit steep dives without exceeding the critical 
Mach Number and to obtain favorable nose-up pitching 
moments, a dive flap was fourtd to be a vital require- 
ment. Investigations in the 8 by 12 ft. Lockheed wind 
tunnel showed that wing dive flaps, while effective, 


’ caused extremely high balancing tail loads. The pres- 


ence of the landing gear, wing flaps, and ailerons limited 
the location of such wing flaps severely. Cowl flap 
arrangements on the aft fuselage caused undesirable 
tail-buffeting and improper trim changes. After 
lengthy testing, the location on the fuselage behind 
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the nose gear was found to be good from all points of 
view. The tail load is reduced with the flap down, the 
drag increment high, and the pitching moment mild 
but favorable, increasing in effect with Mach Number. 
Fig. 33 shows the effect of the flap in increasing the dive 
angle at a Mach Number of 0.77. Without the flap, 
a dive angle of only 8° can be obtained at 30,000 ft. 
with power off before accelerating to high Mach Num- 
bers. The dive flap allows the glide angle to be in- 
creased to 18° for the same speed. 

The speed of operation of the flap (2 sec.) makes it an 
excellent aid in formation flying. Its use on take-off 
helps to lift the nose gear off the ground and benefits 
maximum lift slightly. Mechanically, it provides rapid 
access to hydraulic and electrical equipment on the 
ground without requiring additional removable panels. 


Icing Problems 


While there are numerous cases recorded of P-80A 
aircraft flying through icing conditions, none indicated 
any appreciable formation of ice on the aircraft or en- 
gine intake. This is probably due to a number of causes, 
the most important being the adiabatic temperature 
rise due to speed (Fig. 34). Even at climbing speeds, 
a temperature rise of over 20°F. is available over much 
of the airplane, including the engine inlet screen. The 
plenum chamber obtains full adiabatic rise. At high 
speed, the temperature rise is between 50° and 70°F., 
so little possibility exists for icing then. It is possible, 
based on experience to date, that high-speed jet air- 
craft have an icing problem during only a short period 
on take-off and landing. Additional research should 
be undertaken at once, however, to confirm the above, 
since there may be conditions of supercooled droplets 
being encountered which would permit bad icing in 
climb. Because of its low thrust and drag, a jet air- 
craft is not basically well adapted to handle much ice. 


FLUTTER AND VIBRATION PROBLEMS 


Ailerons 


The most interesting vibration problem encountered 
during the development of the P-80 was that of aileron 
“‘buzz’”’ and oscillation. At Mach Numbers close to 
0.80, the aileron trailing edges would develop a high- 
frequency buzz with an amplitude of about '/,in. At 
higher speeds the whole aileron would oscillate violently 
through large amplitudes. Flight tests were under- 
taken to determine the cause of these phenomena and 
soon established the facts that it was not a case of or- 
dinary flutter and that compressibility effects on the 
wing were at the bottom of the trouble. 

Lockheed then constructed and instrumented a test 
rig consisting of a full-scale outer wing panel and 
aileron, the aileron booster and hydraulic system, the 
control stick, and other equipment items. This ar- 
rangement was then installed in the N.A.C.A. 16-ft. 
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diameter high-speed wind tunnel at Moffet Field and 
tested. An intensive program was undertaken which 
confirmed the flight-test conclusions and furnished 
much other valuable data. The basic cause for the 
aileron movement was shown to be the pressure varia- 
tions and turbulence behind the shock waves that 
formed on the top and bottom of the wing. The motion 
of the shock waves fore and aft on the wing was af- 
fected or aggravated by the movement and deflection 
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IMPORTANT ITEMS TESTED IN HIGH SPEED TUNNEL 
TO PREVENT AILERON OSCILLATIONS 
DESCRIPTION SKETCH COMMENT 
INCREASE ¢ DECREASED 
4 \sraric DYNAMIC BALANCE \ 
INCREASED WING TORSIONAL 
2. | (ADDED A1BS) WO. EFFECT 
LEADING EDGE COUNTER- 
4. |RIGID TIE-DOWN OF TIP Kit NO EFFECT 
SPO“ ERS MANY VARIATIONS OF ANTI 
5. AILERON ¢ WING CONTOUR r alt. SERVO TABS, AILERON SECTIONS 
L VENTS, SPOLERS ETC. 
| PROVED INEFFECTIVE 
SUCCESSFULLY DAMPED OSCTILLATIO- 
6. — NS WITH LOOSE CABLES TO HIGH- 
EST AVANABLE TEST MACH NO(83) 
\WIGH AILERON CONTROL CABLE 350" Tension | successfuL Im ¢ IN 
7. |TENSION UTILIZING STANDARD AiR- SERVICE TO HIGHEST MACH. 
PLANE WYDROULIC CONTROL BOOSTER NUMBERS EXPERIENCED 
Fic. 35. 


Fic. 36. Three types of Lockheed acceleration-type flutter 
dampeners. 


of the aileron itself which occurred under the high air 
loads. Fig. 35 gives a brief summary of the important 
items tested and their effect on the forced oscillations. 
Changes in static and/or dynamic balance were inef- 
fective as were increases in wing torsional stiffness and 
changes in the location of the center of mass of the 
wing. Rigidly restraining the wing in torsion and 
bending was also completely ineffective. Along aero- 
dynamic lines, efforts to fix the location of the shock 
wave by spoilers on both surfaces, venting, changing 
the aileron contour, and making miscellaneous use of 
all types of tabs were also ineffective. Angle-of-attack 
variations of several degrees did not play an important 
part in determining the critical Mach Number for os- 
cillations to begin. The above tests were all run with 
control cables disconnected. 

It was soon noted that aileron cable tension was ex- 
tremely effective in stopping aileron oscillation to the 
highest available test Mach Number. The combination 
of friction and the damping from the hydraulic booster 
provided an effective solution to the problem. Flight 
tests showed that, because of thermal differential con- 
traction of the wing and the control cables, the control 
rigging would decrease from 150 Ibs. to zero at 25,000 


1947 


ft. By the mere expedient of making the aileron ribs 
strong enough to withstand the “‘buzz’’ condition and 
rigging the control cables to 350 Ibs. or more to prevent 
oscillations, a practical solution to the problem was ob- 
tained. No service reports on aileron oscillation have 
been received with the above changes in effect. Cable 
tensions of 50 Ibs. or more were sufficient to stop the 
oscillations, so a wide margin exists witha sea-level rig- 
ging of 350 lbs. Tests were also run to measure the ex- 
citing force, which was found to be small and also, when 
once started, would not increase with Mach Number to 
any important degree. 

Prior to encountering the aileron problem on the P- 
80A, Lockheed had been doing research on flutter 
dampeners. Mr. Bleakney and Mr. Beman were suc- 
cessful in developing an acceleration-type dampener 
(Fig. 36), which, when tested in the wind tunnel, suc- 
ceeded in preventing aileron oscillation with the control 
cables disconnected. It has not been used in flight 
to date because of the simpler solution available for the 
P-80A, but it is being installed on an aircraft for test in 
the belief that such a unit will have definite use on fu- 
ture aircraft. 

As a result of the tests, it is the writer’s opinion that 
airfoil section (particularly position of the maximum 
ordinate), trailing-edge angle, and thickness are the de- 
termining factors as to whether forced aileron oscilla- 
tions will be obtained at the critical Mach Number. 
Some flight-test verification of this has been obtained. 
Also, the oscillating pressures would tend to make the 
use of aerodynamic balance at high speeds question- 
able. 


Elevators 


No flutter problem was encountered on the elevator 
in its normal configuration. One aircraft was lost be- 
cause of failure of a trim-tab mechanism under fatigue 


Fic. 37. Full-scale horizontai tail installed in Lockheed wind 
tunnel for elevator flutter tests with freetab. Strut with bumper 
limits tip deflection. 
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loads. With a free tab, violent oscillations of the eleva- 
tor caused stabilizer failure. Wind-tunnel tests (Fig. 
37) on an actual tail surface with a free tab confirmed 
the speeds at which flutter could be obtained and the 
great violence of this type of flutter. 


Rudders 


No vibration problems have been encountered on the 
rudders. It is, however, extremely difficult to keep rud- 
der counterweights on in spins. Apparently, a 50g 
analysis condition is not sufficient to withstand the 
buffeting in this maneuver, and several balances were 
lost in testing. The pilot was unaware of this loss in 
each case, and no flutter developed in spite of the loss of 
the balance weight. 


ARMAMENT AND CAMERAS 


Armament 


Much of the P-80 development has been concerned 
with the various types of armament installed and in 
studies of means of protection of the engine installation 
from battle damage. Because of the close concentra- 
tion of six powerful guns in the nose, a considerable 
structural problem presented itself in making the for- 
ward nose of the fuselage and the gun blast tubes take 
the beating that occurred. The problem was finally 
solved after extensive ground and air firing, in which 
several aircraft shot more than three times the number 
of rounds averaged by fighter aircraft in the war. The 
disposal of empty cases and links to prevent their en- 
tering the ducts was solved by a retractable chute and 
proper detail design and shielding of critical items in 
the plenum chamber. The entrance of corrosive powder 
gases into the ducts and engines requires special pre- 
cautions by painting or servicing in the field. This is a 
common problem on several fighter types. 

Fig. 38 shows the aircraft equipped with two 1,000- 
Ib. bombs. With the exception of the previously men- 
tioned problem of wing-tip vortexes, no special problems 
have developed with this installation. 

An interesting design of rocket launchers on the P-80 
is shown on Fig. 39. 

A total of ten large rockets may be carried under the 
wings on a special rack designed for installation or re- 


Fic. 38. P-80A carrying two 1,000-lb. bombs. 
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Cluster of five rockets (without warheads) hung 
under each wing of P-80A. 


Fic. 39a. 


Fic. 39b. Nose camera installation showing nose hinged upward. 


moval in 30 min. The aircraft is capable of doing all the 
various tasks expected of fighter types, such as strafing, 
bombing, ground attack with rockets, interception, 
bomber escorting, and photographic work. 


Photographic Version 


Fig. 39a shows the installation of a special fuselage 
nose containing various cameras. Standard fighters 
can be readily converted into photographic ships by 
interchanging this section of-the aircraft. The ready 
accessibility to the cameras provided by hinging the 
nose forward is apparent. 


Battle Damage Tests 


The vulnerability of turbojet engines in fighter air- 
craft was unknown in 1944. The fuselage installation of 
the I-40 engine on the P-80A presented a number of 
special problems regarding liberated heat with a punc- 
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Fic. 40. Official A.A.F. photo of YP-80A in Wright Field wind 
tunnel arranged for battle-damage tests. 


Fic. 41. Official A.A.F. photo of wind-tunnel test set up to 
shoot YP-80A under operating conditions in Wright Field 20-ft. 
wind tunnel. Note 0.50-cal. machine gun to right of armor plate 
used to stop damaged turbine wheel or other engine parts. 


tured tailpipe or flame tubes. To obtain answers to 
these problems, Lockheed proposed to the A.A.F. that 
a series of investigations be undertaken in a wind tunnel 
involving actual shooting of an operating engine-air- 
plane combination. This program was accepted and 
carried out by the A.A.F. in the 20-ft. high-speed wind 
tunnel at Wright Field under the supervision of Col. 
George Price. 

Figs. 40 and 41 show the test setup used. Great haz- 
ard to the tunnel existed during the test, since the engine 
was shot in various places when operating at imilitary 
power and at actual airplane flight speeds. Invaluable 
information was obtained leading to greatly reduced air- 
craft vulnerability. Results of these basic tests will be 
of great assistance in proper design of later military air- 
craft. 
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PROBLEMS OF PILOT COMFORT 


Accelerations 


Because of its high speed and maneuverability, the 
P-80 presents a considerable problem of pilot comfort 
at high accelerations. In level flight at low altitude on 
a hot day, accelerations of 4g to 5g are constantly being 
indicated by the accelerometer. These are sharp, 
solid bumps and are uncomfortable. The safety belt 
must be kept tight, and a crash helmet to prevent 
bumping the head on the canopy is definitely required. 
The use of g-suits is necessary to obtain the maximum 
maneuverability out of the airplane, but the net result 
is merely to raise the operating limits of the pilot closer 
to those of the aircraft. It is distinctly possible that the 
use of g-suits will require higher design fighter load fac- 
tors. 


Cockpit Temperature Control 


The effect of high speed on cockpit air temperature 
due to adiabatic rise is shown on Fig. 34. Solar radiation 
through the canopy adds about 30°F. to the air tem- 
perature, making high-speed flight at low altitude 
extremely uncomfortable. It is discouraging to turn on 
some nice fresh air to cool your face and find out that a 
hot blast comes out of the air duct. 

Pressurization of the cockpit using engine compres- 
sor air increases the cooling requirements as shown on 
Fig. 42. The cockpit pressure is kept to a moderate 
amount because of rapid decompression that would re- 
sult from battle damage to the canopy at high altitude. 
Cockpit pressure is turned on at 8,000 ft., resulting in 
the pressure curve shown and a rapid temperature in- 
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crease of an unbearable amount unless cooling is used. 
Line A-B is the cockpit air temperature in climb (ac- 
counting for the adiabatic rise in climb and a small in- 
crease for the solar effects at the time noted) before 
pressure is turned on. At 8,000 ft., the cockpit air tem- 
peratures would increase along line B-C and then de- 
crease along C-D with altitude and no cooling. It is 
practical to pressurize only above 20,000 ft. in this 
case. 

To obtain proper pilot comfort, the unit shown in Fig. 
4 was developed to cool the cockpit air. An inter- 
cooler, turbine, and fan decrease the inlet cockpit air 
as shown by the ‘Controlled Temperature Range.”’ 
This cooling unit has been extensively tested in the P- 
80A and is now in production by AiResearch. 


Safety Features 


Machmeters are definitely required for safe operation 
of jet fighters. This instrument has proved reliable in 
service and is a standard item of equipment. 

Fig. 43 shows a mock-up of a jettisonable seat de- 
veloped along the German lines to provide better es- 
cape means at high speed. The close clearances for in- 
stalling such a seat are readily apparent in the photo- 
graph. This installation will be available in service 
on the P-80B. 


FLIGHT TESTING 


Equipment 


In general, flight testing of jet aircraft tends to be 
simpler than that of conventional types. The simplicity 


Fic. 43. Mock-up of jettison seat installation. 
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Fic. 46. Nose support and adjustment on ground thrust stand. 


of the engine installation, as well as the lack of power 
effects on stability, lift, and drag making application of 
wind-tunnel data fairly easy and accurate, tends to 
eliminate much testing required previously. The type 
of testing becomes more along aerodynamic lines in 
that pressure surveys, temperature distributions, and 
velocity measurements become more important in the 
test program than the development testing of a maze of 
mechanical gadgets. 
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Fic. 48. P-80A taking off assisted by two 1,000-lb. thrust 
rocket units. 


The most useful item of equipment developed during 
the P-80 testing is shown in Figs. 44, 45, and 46. It is 
a thrust stand that allows the complete airplane to be 
suspended on a cable system and the engine thrust 
measured accurately on a scale. Because no cooling 
problem of an appreciable magnitude exists during 
ground running, a substantial amount of installation 
and engine development can be undertaken with simple 
equipment. The accuracy of the equipment shown is 
good, being of the order of '/2 of 1 per cent. It cannot 
be used in the wind, but in several cases the aircraft 
has been set up in a hangar and run with the tailpipe 
pointed out a door, so this is not much of a limitation. 
In several tests, engines calibrated at the factory and in 
the thrust rig have given data on thrusts and fuel flows 
checking within 1 per cent. During all accurate per- 
formance tests, thrust measurements are taken before 
and after the flights to be sure that engine performance 
has not deteriorated. Because of the great variation of 
performance between production engines, as well as a 
decrease of thrust approaching 14 Ibs. per hour of oper- 
ation, it is essential to do this to get satisfactory testing 
accuracy. 

Another useful arrangement developed for P-80 
testing is a two-place modification to the cockpit allow- 
ing a flight-test engineer to accompany the pilot to take 
additional data and observe the tests (Fig. 47). It is 


believed that about 30 per cent of the flight time can be 


saved using this installation, in addition to giving the 
engineers first-hand knowledge of the problems en- 
countered. 

Fig. 48 shows the application of rocket assist for take- 
off at heavy loads. While this means is effective for re- 
ducing take-off distance, the cost of $250 for a 30 per 
cent reduction cannot be considered justified for normal 
operation. Water injection reduces the take-off dis- 
tance 20 per cent much more economically. 

Thrust measurements in flight involve the use of a 
tailpipe survey rake shown in Fig. 49. This instrument 
is useful and gives good results when combined with 
other equipment measuring engine inlet conditions. 

Air-speed measurements, especially in dives, require 
careful balancing of the total pressure and static lines, 
as well as small volumes in these systems. In spite of 
taking the greatest pains with the flush static installa- 
tion on the fuselage nose, the air-speed calibration 
varies a considerable amount from airplane to airplane 
(sometimes 1 to 2 per cent). This requires special 
calibration of each test aircraft, usually by the altimeter 
method at low altitude. 

Fig. 50 shows a typical P-80 test instrumentation in- 
stallation for automatic data recording. 


Radio Control Drones 


It is increasingly evident that radio-controlled drone 
planes will be invaluable for high Mach Number re- 
search, testing jettisonable seats, and structural inves- 
tigations to high load factors or in sudden maneuvers. 
A number of P-80A airplanes incorporating equipment 
developed by Bell, the U.S.A.F, Sperry, and others have 
been built. Experience soon should be obtained with 
their use. 


Fic. 49. Tailpipe survey rake used to determine thrust in flight. 
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PROBLEMS IN SERVICE 


Power Plant 


It is natural that service experience with such a new 
type as the P-80, involving so many new features, 
would bring to light a host of new problems. The lack 
of propeller slipstream on take-off and in low-speed 
maneuvers has been confusing to pilots trained in con- 
ventional aircraft. The lack of dependable means for 
preventing engine overspeeding has overstressed tur- 
bine wheels in many cases. Even in cases where over- 
speeding has not been present, turbine wheels have 
failed with disastrous results. Figs. 51 and 52 show 
what happened in one case where a turbine wheel 
failed in flight close to the ground. In this case, the 
pilot landed safely after the wheel split and flew through 
the sides of the fuselage. During an early Lockheed 
test, Tony LeVier had a wheel fail at 540 m.p.h. This 
cut the whole aft end of the airplane off. While the 
pilot escaped safely, the turbine wheel was not found 
for a year, even after a most extensive and highly 
scientific search. A farmer finally picked up a large 
piece 6 miles to one side of the flight path. There is no 
practical means of providing a guard to absorb the 
energy from a large split turbine wheel to prevent it 
from doing extremely great damage. The wheel must 
be designed to be practically infallible even after losing 
a number of blades. So far, blade losses have been few 
and not at all dangerous. They generally go out the 


Fic. 51. Crash landing after turbine-wheel failure showing 
damage to left side of fuselage. 


Fic. 52. Damage to right side of fuselage due to turbine-wheel 
leaving in flight. 
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Fic. 53. Damage to paint on fuselage and fillet occurring in 11.3g 
pull-out at high speed. 


Fic. 54. Effect of erosion on nose and wooden dummy guns after 
flying at high speed in rain for 2 min. 


tailpipe. Reliable governing is an absolute must to 
prevent wheel failures. 


High Accelerations 


The high maneuverability and speed of the P-80 
and the use of g-suits has resulted in the aircraft being 
subjected to a severe beating in service. Fortunately, 
this problem was anticipated to some degree, and in 
static tests the wing took 30 per cent over its design 
load factor. The other parts of the airplane, particu- 
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Fic. 55. Effect of hitting rain for 2 min. at high speed—wing 
leading edge. 


larly the tail and fuselage, were also overstrength. The 
highest load factor recorded (on an N.A.C.A. V-G 
recorder) on a P-80 has been 11.3g, which occurred dur- 
ing gunnery tests too close to the ground. 

Fig. 53 shows the damage to the fuselage fillet and 
paint cracking at joints. Wing leading edges had to be 
replaced. In tests, over 9g acceleration has been held 
throughout a 360° turn. The occurrence of 8g in service 
has been common. 

With such load factors being obtained, the only 
structural trouble encountered has been in nonstressed 
(?) fillets and trailing edges where flush riveting has 
been used. The use of paint on the surface, however, 
creates an unfavorable impression, particularly when it 
cracks at seams and over rivets. 


Aircraft Finish 


This subject has been covered somewhat, but it 
might be added that maintenance of aircraft finish in 
the field is extremely difficult. Between the load fac- 
tors obtained by exhuberant pilots, hob-nailed-booted 
mechanics, desert sandstorms, rain, and general service 
use, the paint job does not have much of a chance. 
Soon the paint is so rough that the airplane is slower 
than if it were not painted. Leading edges roughen 
up and the stall characteristics deteriorate, most gener- 
ally with wing-tip fuel tanks on. 

Figs. 54 and 55 show paint and wood (dummy gun) 
erosion after 2 min. in a heavy rain. 

The only immediate solution to this problem is to re- 
move the paint, but for the gains involved further 
finish development is certainly in order. 


Fueling 


Because of the nature of its fuel system, the fuel for 
jet engines must be filtered with micronic filters. All 
fuel trucks servicing jet aircraft must be so equipped 
(as must the engine itself), and pains must be taken to 
be sure no water is mixed with the fuel. At high 
I.p.m., a jet engine can swallow a good deal of water 
in the fuel without danger, but in landing with low 
r.p.m. it is dangerous because it can put the flame 
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LOCKHEED P-80A JET FIGHTER AIRPLANE 


Fic. 56. Refueling and servicing a P-80A in 4 min. 


Fig. 56 shows a record fueling-up time being estab- 
lished for a P-80. It took slightly over 2 min. to put 
in 750 gal. of fuel. 


Navigation 


With jet engines, the fuel consumption at sea level is 
about three times that at 30,000 ft. at a given r.p.m. 
This results in short sea-level endurance which, com- 
bined with a cruising speed of well over 8 miles per 
min., makes it hazardous and easy to get lost. The pilot 
must do an excellent job of navigation under these con- 
ditions. In bad weather when low-altitude flying is re- 
quired, the problem is greatly accentuated. Better 
radio facilities, on the ground as well as in the aircraft, 
must be obtained. The losses of aircraft due to the 
above causes inevitably run higher than with conven- 
tional aircraft. 


CONCLUDING COMMENTS 


The experiences encountered during the development 
of the Lockheed P-80 “‘Shooting Star’ fighter lead the 
author to attempt to summarize what he believes to be 
the outstanding problems requiring solutions to improve 
the safety and utility of jet-powered aircraft. The aero- 
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Mayor Jer AircRAFT DEVELOPMENT PROBLEMS 


|. BETTER NAVIGATION AND LOW W/S/BILITY LANDING 
FACILITIES. 


2. /MPROVED TAKE-OFF PERFORMANCE FROM BETTER 
WINGS, FLAPS, OR BOUNDARY LAYER CONTROL. 


3. /MPROVED CREW ESCAPE FACILITIES AT 
HIGH SPEEDS. 


4. IMPROVED LANDING GEARS, TIRES, ANO BRAKES 
(OR OTHER MEANS FOR DECELERATION). 


5. RELIABLE ANTI-FLUTTER DEVICES. 


6. RELIABLE FIRE DETECTION AND EXTINGUISHING MEANS. 


7. ELECTRIC WINDSHIELD AND CANOPY DEFROSTERS 
AND DEICERS. 
8. SUITABLE FILLERS, PUTTIES, AND PAINTS TO 


OBTAIN AERODYNAM/C SMOOTHNESS. 


Fic. 57. 


Mayor Jer Encine Devetopment Proe.ems 


BETTER OVER-ALL RELIABILITY. 


READ DUAL FUEL SYSTEMS. 


MORE TAKE-OFF THRUST. 


4 
3. LOWER FUEL CONSUMPTION. 
4. 


BETTER ACCELERATION FROM LOW RPM.» 
Now AT HicH &.2M.) 


MORE RELIABLE GOVERNORS. 


N 


BETTER POWER CONTROLS. 


8. A ORY AIR PRESSURE SOURCE WHEN 
USING FLUID INJECTION. 


9. LESS FUEL LEAKAGE FROM THE ENGINE 
ITSELF ON STARTING. 


10. BETTER STARTING CHARACTERISTICS 
PARTICULARILY AT ALTITUDE. 


Fic. 58. 


dynamic and power-plant improvements required to 
get higher performance are well known and are not in- 
cluded in the tables. An attempt has been made to put 
the various items in the order of their importance. 
Figs. 57 and 58 are self-explanatory. 

There seem to be enough problems left to keep us all 
busy for some time to come. 
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Approximations Involved in the Linear 
Differential Equation for Compressible Flow 


H. W. SIBERT* 
Unwersily of Colorado 


SUMMARY 


An investigation is made of the approximations necessary to 
linearize the partial differential equation defining the velocity 
potential for compressible flow. It is found that not only second- 
order terms but also a first-order term must be negligibly small 
compared with unity. The application of these conditions to the 
Ackeret theory for two-dimensional supersonic flow gives both 
upper and lower limits to the free-stream Mach Number for 
which this theory will give reliable results. 


SyMBOLS 


In the discussion that follows it will be assumed that the free- 

stream velocity is in the direction of the positive x axis. Let 
u,v = components of velocity in direction of x, y axes 

= velocity potential (u = ¢:,0 = ¢,)f 

speed at any point = V/ u? + v? 

free-stream speed 

u—U 

speed of sound at any point 

ratio of specific heats 

static pressure at any point 

dynamic pressure at any point 

(pb — po)/qo = pressure coefficient at any point 

0 = U/ao = free-stream Mach Number 

Norte: Subscript 0 refers to free-stream conditions. 


INTRODUCTION 


Fe THE SAKE OF BREVITY the following discussion 
will be limited to a two-dimensional flow in planes 
parallel to the xy plane. The extension of the discus- 
sion to the three-dimensional case can easily be made. 

When the flow is steady, isentropic and irrotational, 
the differential equation defining a two-dimensional 
flow in planes parallel to the xy plane is 


(1 — — v?/a*) by — 2(uv/a*)dry = 0 (1) 
and the pressure coefficient at any point of the flow is 


P = {1 + — — V2/U2) P 
Q) 


It is often stated that, when second-order terms are 
neglected, Eqs. (1) and (2) become, respectively, 


(1 M_’)¢,. + ¢y = 0 (3) 
P = —2u’'/U (4) 
Apparently, no one has yet defined exactly what ap- 


Received March 6, 1947. 

* Professor of Aeronautical Engineering. 

¢ In the subscript form of partial derivatives ¢. = 0/dx, 
gry = etc. 


proximations must be made to obtain Eqs. (3) and (4) 
from Eqs. (1) and (2). The determination of these 
approximations is the purpose of this paper. 


DETERMINATION OF APPROXIMATIONS 


The most obvious assumption to make is that the 
velocity at any point does not differ much from the 
free-stream velocity. This can be accomplished by 
assuming that: 

(1) (u’/U)? and (v/U)? are negligibly small com- 
pared with unity. 

Since u = U+ u’ and V? = u? + v?, it follows from 
this assumption that, to a close approximation, 


= ut/U? = 1 + u'/U) (6) 

Then Eq. (2) becomes 
P = (2/yM.2){[1 — (y — I) 
(6) 


By means of a binomial expansion Eq. (6) will reduce 
to Eq. (4) if the following assumption is made: 

(2) [(y — 1)Mo?(u’/U)]? is negligibly small com- 
pared with unity. 

From the Bernoulli equation, 


a? + (y — 1)(V?/2) = a? + (y — 1)(U?/2) (7) 


By means of Eq. (5) and the relation Mp = U/ao, Eq. 
(7) becomes: 


= 1 — (y — (8) 
which by virtue of the second assumption reduces to 
a*/ao? = 1 (approx.) (9) 


Then, from Eqs. (9) and (5) 
U?/a* = (U/ao)?(ao/a)? = Mo? (approx.) (10) 
u*/a? = (u/U)*(U/a)*? = [1 + 2(u’/U)]Mo? (approx.) 
(11) 
v?/a® = (v/U)*(U/a)? = (approx.) (12) 
uv/a? = (uv/U*)(U/a)? = M,*uv/U? (approx.) (13) 
It follows from Eqs. (11) and (12) that the first two 
terms of Eq. (1) will reduce to the left side of Eq. (3) 
if it is assumed that: 
(3) 2u’/U and (Mw/U)? are negligibly small com- 
pared with unity. From this assumption, the relation 
u = U+u’' and Eq. (13) 
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EQUATION FOR COMPRESSIBLE FLOW 681 


u/U = 1+ (u’/U) = 1 (approx.) (14) 
uv/a? = My*v/U (approx.) (15) 


Hence, in order for the third term of Eq. (1) to be 
neglected, it must be assumed that: 

(4) 2Mo?(v/U)¢z, is negligibly small compared with 
either (1 — or 

Since (u’/U)? < 2u’/U when the latter is small com- 
pared with unity, the above four assumptions can be 
written as follows: 

(a) 2u’/U, (v/U)*, (Mw/U)? and [(y — 1)M? X 
(u’/U)]? are negligibly small compared with unity. 

(b) 2M,?(v/U)dzy is negligibly small compared with 
either (1 — Mo?) OF dyy- 

It is evident from the first term of assumption (a) 
that a first-order term, in addition to second-order 
terms, must be neglected in order to obtain Eqs. (3) 
and (4) from Eqs. (1) and (2), respectively. 


APPLICATION TO ACKERET SUPERSONIC THEORY 
By introducing the disturbance potential ¢ defined 
by the relation 
¢=¢+ Ux (16) 
Eq. (3) becomes 
(1 — + = 0 (17) 
and the disturbance velocity components u’ and v be- 
come 
(18) 
It can be shown that, when the positive y axis for the 
upper surface of the airfoil points upward and the 
positive y axis for the lower surface points downward, 
the correct solution of Eq. (17) is: 
r=x-VMi- 10) (19) 
where f is any function of r. Note that, from Egs. 
(18) and (19), 


v/u! = ¢,/¢2 = —1 


ul = $2, v= 


(20) 


If the airflow follows the contour of the airfoil, the 
slope y’ (or dy/dx) of the airfoil becomes [see Eg. (14) ]: 


y’ = 0/u = (v/U)(U/u) = (approx.) (21) 
It follows from Eqs. (20) and (21) that 
u’/U = (u'/v)(v/U) = —y’/V Mo? — 1 (approx.) (22) 


By virtue of Eq. (21) and (22), assumption (a) of the 
last section becomes 

(A) 2y'/V Me? — 1, (9')*, (May’)?, and [(y 1) X 
M,*y'/V M,*? — 1]? are negligibly small compared with 
unity. 

At supersonic speeds the second term of assumption 
(A) is unnecessary because it is less than the third 
term. Since y = 1.4 approximately, the fourth term of 
assumption (A) can also be neglected because it is less 
than the first term when My) < 2.2 and less than the 
third term when My, > 1.1. Thus, assumption (A) 
reduces to: 

(A’) 2y’/ VM? — 1 and (Moy’)? are negligibly small 
compared with unity. It follows from assumption 
(A’) that 

(1) Everywhere on the airfoil, y’ must be small com- 
pared with unity, which requires a thin airfoil with 
sharp leading and trailing edges at a small angle of at- 
tack. 

(2) Mo must not be close to unity. 

(3) (Moy’) must be small compared with unity. 


Hence, the higher the Mach Number, the thinner 
the airfoil and the smaller the angle of attack if accurate 
results are to be obtained from the Ackeret theory. 
Thus, for any airfoil at a given angle of attack, there is 
both a lower and an upper limit of the value of M) for 
which Ackeret’s theory will give reliable results. 

The validity of the Ackeret theory below My = 1.4 
and above My) = 4.5 is questionable. Hence, this 
theory should be used with caution outside of this 
range until sufficient test data are available to indicate 
what range of Mo should be used for any particular 
airfoil section at a given angle of attack. 


Letter to 


Dear Sir: 

The excellent paper on the efficiency of gas-turbine power 
plants for stationary and jet propulsion purposes by Prof. J. H. 
Keenan and Prof. J. Kaye in the August, 1947, issue of the 
JourNAL is certainly much welcomed by the technical’ circles 
concerning the design of such power plants. 

The value of the curves can be extended somewhat by more 
suitably choosing the parameters. The writer made some similar 
investigations on turbojet and ram-jet in 1945 and derived, in an 
unpublished paper, a general formula of the propulsion effi- 
ciency, under the assumption of ideal gas with constant specific 
heat and negligible fuel-to-air mass ratio. The formula is: 


the Editor 


tT — ty Ir — &@ — 0) — (&/r)(e — 
_vizi} 


—Vé-— 1 is the ratio of the flying speed to the maximum gas 
speed obtained by expanding atmospheric air to vacuum; this 
ratio is the Mach Number of flying speed times VJ/ (k — 1)/2, 
k being the isentropic exponent. + is the ratio of maximum 
temperature at turbine entrance, or in the case of ram-jet, at 
entrance to jet nozzle, to atmospheric temperature. ¢ signifies 
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the compression ratio of the compressor. In case of turbojet, 
it is defiged as the ratio of exit stagnation temperature to suc- 
tion stagnation temperature of compressor, and it is reduced to 
unity in case of ram-jet. The losses are taken care of in the 
quantities: 


= n2 = = — 1 


The n’s are the efficiency of different processes, the indexes r, c, 
t, n indicate the ram effect, compressor, turbine, and jet nozzle. 
np takes care of the pressure losses and is defined as (P,;P.)%» = 
(P.P,), where the P’s are the pressure ratios (>1) and the indexes 
have the same meaning as above. All efficiencies are polytropic 
efficiencies—i.e., relative to the theoretical work along the actual 
line of change of states, instead of isentropic. For compression 
shock at supersonic speed 1 — 7, is the ratio of actual increase 
of entropy to isobar increase of entropy between the same 
temperatures. 

The advantages of such representation are evident: 

(1) The number of independent variables is reduced to a 
minimum. The actual atmospheric conditions and flying speed 
are immaterial because only ratios appear in the formula. 

(2) By using polytropic instead of isentropic efficiencies, the 
variation of machine performances with compression ratio is 
avoided. 

(3) Analytical treatment of the relatively simple formula is 
possible. For example, it is quite easy to find the optimal tem- 
perature (r) for the ram-jet by putting the derivative of » with 
respect to r to zero. It will be found that, other things being 
equal, the optimum temperature depends on the speed, and there 
is no reason to use very high temperature even if the material 
can stand it, In case of turbojet, high temperature is, of course, 
desirable because one can always adjust the compression ratio 
to suit. 

For the stationary power plant similar results can be obtained. 
Taking the general case of n-stage compression and m-stage ex- 
pansion with regeneration and assuming recooling and reheating, 
respectively, to initial stagnation temperature of compressor 
intake and turbine inlet, we obtain the thermal efficiency of the 
cycles: 


mr(1 — — — 1) 


or(1 — — — 1) 
or(1 — E-7/%) + — 


where ¢ = m/n. 1 is the ratio of stagnation temperature at 
turbine inlet to compressor intake. ¢ signifies the compression 
ratio and is defined as g = £, where £ is the ratio of stagnation 
temperature at compressor exit to compressor inlet as above. 
n’'=nenmp.* (1—¢) corresponds to the regeneration efficiency as 
defined by Keenan and Kaye and {’ = {/n. The advantages 
over the conventional formula are the same as for jet propulsion, 
and one notes especially the small number of independent 
variables. 

The effect of variable specific heat can also be taken care of 
with one single additional parameter. Defining Cp and Cp’, the 
linear and logarithmic average of specific heat, respectively, 
between compressor inlet and outlet temperature—i.e., 


— Ti) = = — 
cp’ log, (T2/T:) = (cp/T)4T = 2 — 


cp and c,’’ for the turbine temperature interval, cp* the linear 
average specific heat between the temperature at outlet of heat 
exchanger and the temperature at turbine exit, it can be shown 
that the above formula becomes 


where ¢ = mCp/Nncp, = = f'cp*/cp. The enthalpy 
h and the entropy ¢ are listed in reference 1 of the paper. 

Even the problem with jacket-cooled compression and ‘‘wet 
compression’”’ (Kleinschmidt, Mechanical Engineering, Vol. 69, 
pp. 115-116, 1947) can be attacked similarly by defining the 
quantities ¢ and o’ in a suitable manner and by assuming steam 
as ideal gas. 

* np is defined by (P-)?? = (P+) as before. 


L. S. DzuncG 
Brown, Boveri & Company Limited 
Baden, Switzerland 
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Free Vibrations of Sweptback Wing 


ADAM T. ZAHORSKI* 
Northrop Aircraft, Inc. 


SUMMARY 


The frequencies and modes of free vibration, coupled in bend- 
ing and twisting, of a free-free, sweptback, uniform beam are ob- 
tained from its differential equations of motion. The modes of 
vibration of a nonuniform, sweptback wing are then determined 
by the Ritz method. The discussion of the methods and various 
practical aspects of the analysis are presented from an engineer- 
ing, rather than a mathematical, point of view. 


NOTATION 


reference rectangular coordinates 

coordinate along sweepback span 

angle of sweepback 

semispan along the sweepback 

flexural rigidity, lb.-in.? 

torsional rigidity, Ibs.-in.? 

distributed mass of wing, Ib.-sec.? per in.” 

distributed mass polar moment of inertia in Ib.- 
in.?-sec.? per in,? 

radius of gyration, mro? = Ip 

m/EI 

GJ/EI at the origin 

V EI/GI(10/b) 

downward deflection of the beam 

angle of twist 

time, sec. 

frequency of vibration, rad. per sec., uniform beam 

frequency of vibration, rad. per sec., nonuniform 


beam 
kW 
bending mode, uniform beam 
twisting mode, uniform beam 
bending mode, nonuniform beam 
twisting mode, nonuniform beam 


normalized modes 
normalizing factor 


on 


~ 
wou 


INTRODUCTION 


—_ THE END OF THE WAR the sweptback wing has 
become popular with most aircraft designers. It 
has also introduced many new aerodynamical and struc- 
tural problems demanding the immediate attention of 
engineers. One of these problems pertains to the 
modes of vibration of the swept wing. In this paper a 
theory is developed for the modes of vibration of a non- 
uniform, sweptback beam, and it is shown how this 
theory can be applied to an airplane such as the XB-35 
Flying Wing (Fig. 1). 

It is generally recognized that in a wing with large 
sweepback the interaction between bending and twist- 
ing may greatly affect the modes and the frequencies 
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of vibration; the problem is complicated and requires 
careful re-examination of the old assumptions and sub- 
stantiation of the new ones. 

The art of mathematical analysis has not as yet pro- 
gressed far enough so that the vibrations of a non- 
homogeneous, three-dimensional body, such as an air- 
plane wing, can be determined. The closest approxi- 
mation to this, for which the equations of motion can 
be set up, is a two-dimensional, nonuniform plate. 
However, our knowledge of the elasticity of plates is 
still limited, and it is doubtful that the necessary elas- 
tic constants for the “‘plate’’ analysis can be caléulated. 

A bent beam, fundamentally one-dimensional “‘line 
beam,”’ possessing properties of mass and stiffness 
(m, I,, EI, and GJ) which vary along its length is the 
best at present that can be assumed for an idealized 
swept wing and for which mathematical analysis is 
plausible (Fig. 2). 

The minimum energy method, known as the Ritz 
method (reference 1, p. 259), will be employed in the 
solution of the problem. This method requires that 
the approximate shape of the modes of vibration be 
known. Any reasonable deformation curves will be 
acceptable provided they satisfy the same boundary 


Fic. 1. XB-35 Flying Wing. 


Fic. 2. Idealized sweptback wing. 
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conditions that exist in the vibrating beam. The 
modes of vibration of a uniform beam with a sweepback 
will be assumed as the reasonable curves on the basis 
of which the Ritz method will be used to obtain the true 
modes and frequencies of vibration of the nonuniform 
beam. 


I.—UNIFORM BEAM 


Consider a uniform, sweptback beam as shown in 
plan view on Fig. 2, where X and Y are reference rec- 
tangular coordinate axes, s is the coordinate along the 
span, and a is the angle of sweep. Assume that the 
distributed mass (m) in Ib.-sec.? per in.? is concentrated 
along the S-axis and possesses the mass polar moment 
of inertia property I, of Ib.-in.?-sec.? per in.?._ Then the 
differential equations of motion of any point (s) are 
given by (reference 1, pp. 222 and 216): 


(0*v/Os*) + = (1) 
(0°8/0s?) — 1°(0?B/ot) = 0 


where 
k* = m/EI and = I,/GJ (2) 


v is the downward deflection, and 6 is the angle of 
twist at the point s. Both v and # are unknown func- 
tions of the time coordinate (¢) and the space coordi- 
nate (s). 

For normal modes of vibration the deflections at any 
point vary harmonically and can be represented by a 
product such as 


v = exp(id)P(s), B = exp(2A/)Q(s) (3) 


where P(s) and Q(s) are functions of the s-coordinate 
only and determine the shape of the normal modes of 
vibration and where ) is the frequency in rad. per sec. 
It must be noted that a mode is represented here by 
two quantities, P and Q, where the P component is 
contributed by the bending of the beam and the Q 
component is contributed by the twisting. For the 
sake of simplicity, P will be referred to as the “bending 
mode” and Q as the “twisting mode.”’ 

If Eqs. (3) are substituted into Eqs. (1), then the 
partial differential equations of deflections v and 8 in 
terms of variables s and ¢ are changed into the total 
differential equations of modes P and Q as functions of 
the variable s, so that 


(d*P/ds*) — = (4) 
(@°Q/ds*) + = 
The general solutions of the above equation are (ref- 
erence 1, p. 223): 
P(s) = Ci cosh ky/Xs + sinh ky/Xs + 
C3cos + Cy sin kx/Xs> (5) 
Q(s) = Bi cos As + Bz sin DAs 


where the six constants must be determined from the 
end conditions of the beam. At the free end of the 


| 
=) 


beam these conditions, which are zero torque, moment, 
and shear, can be expressed as follows: at s = }, 


dQ/ds = d°P/ds? = d’P/ds* = 0 (6) 


Additional conditions can easily be formulated by 
taking advantage of the symmetry of the structure. 
It can be shown’ that in a symmetric structure the 
natural modes of vibration will be symmetric or anti- 
symmetric with respect to the line of symmetry. 

If the beam is cut along the line of symmetry (X-axis 
on Fig. 2) and only the right half is considered, then 
the beam will be in dynamic equilibrium if the section 
at the cut is subjected to the same restraints as in the 
original condition. The two cases are defined as fol- 
lows: 

In the symmetric case only an external moment, My 
about the X-axis (Fig. 3), is necessary for the equilib- 
rium. In this case M) cos a is the moment component 
that resists bending of the beam and —M) sin a is the 
moment component that opposes twisting of the beam 


Me_ Mo cos 
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Fic. 3. Symmetric case. 


To SIN & 


Fic. 4. Antisymmetric case. 


In the antisymmetric case only a vertical shear and 
an external moment 7) about the Y-axis (Fig. 4) are 
necessary for the equilibrium. In this case Ty sin a is 
the moment component that resists bending of the 
beam and 7> cos a is the moment component that op- 
poses twisting of the beam. The additional boundary 
conditions necessary for the solution of Eq. (5) will be 
discussed separately for symmetric and antisymmetric 
cases. 


Symmetric Case 


The conditions at the origin are that there is no ver- 
tical shear, no component of moment reactions in the 
Y-direction, and no rotation of the tangent parallel to 
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the Y-axis. The last two conditions can be easily recog- 
nized on Figs. 5 and 6, and all three can be expressed as 
follows: ats = 0, 


@P/ds* = 0 
EI (d?P/ds*) sin a + GJ(dQ/ds) cosa =0> (7) 
(dP/ds) cosa — Qsina = 0 


Applying the boundary conditions of Eqs. (6) and 
(7) to Eqs. (5), the following relations between con- 
stants is obtained: 


B, = cot a 
By, = B, tan 


ll 


CG, ‘aie C3 = — Bono(l/k?) cot a (8) 


cosh sin + sinh cos (9) 
cosh ~ cos § — sinh £ sin § — 1 


cosh £ sin + sinh cos 


where 
and m=GJ/EI (10) 


at the origin. 

There are six equations from which six constants at 
integration should be determined, but it can be shown 
by successive substitution of constants that in the final 
result all of them will be eliminated leaving only an 
equality 

1 — cosh cos cot A 


—no cot? (11) 


The beam can vibrate only if the above equation is 
satisfied. The solution of this equation, known as the 
“frequency equation” yields an infinite number of 
frequencies \’s, each corresponding to a particular 
mode of vibration. 

It is of interest to note two limiting cases of Eq. (11): 
one when a = 0 and another when a = 7/2. 

In the case of a beam with no sweepback, a = 0 and 
Eq. (11) will be satisfied: (a) if tan Ab = O or lAd = =, 
= 27, = 32, etc. (these are uncoupled torsional 
modes); (b) if 


cosh ky/Xb sin + sinh ky/Xd cos = 0 
or 
tan kx/Xb + tanh = 0 (12) 


This corresponds to the frequency equation for the sym- 
metric vibration of free-free beam given in reference 1, 
p. 277. 

In the case when a = 2/2, Eq. (11) will be satisfied: 
(a) if cot Ab = O or Ab = x/2, = 3/27, etc. (these are 
the uncoupled torsional modes); (b) if 


Fic. 5. 


Fic. 6. 


1 — cosh ky/Xd cos ky/Xd = 0 (13) 


This tallies with the frequency equation for both sym- 
metric and antisymmetric modes of vibration of the 
free-free beam given in reference 1, p. 232. 

It is of interest to observe here that if a beam B’AB 
in X-Y plane (Fig. 7) vibrates in the symmetric modes, 
then its projections on Y-Z and X-Z planes will vibrate 
in such a manner that they will satisfy Eqs. (12) and 
(13), respectively. 

We have seen that all six constants in Eqs. (5) could 
not be obtained. This should be obvious if we consider 
that the amplitude of the modes depends on the amount 
of energy involved in vibration and that this amount 
has not been specified in the problem. On the other 
hand, five of the six constants can be expressed in terms 


Fic. 7. Projections of symmetric vibration. 


is 
685 
pment, 
-2 
b, q 
cos a 
(6) |e 
A G Y 
r anti- ds ] 
Y 
| 2p 
LE i 
, then 
ection 
in the x ee 
is fol- — 
uilib- 2 aq 
onent 
is the | 
ds seco Y 
eam —_ 
"SRK 
ss 7 & 
-Ogs cot a Gy $2 
cos a 
x 
4 
and 
are 
a is 
the A y 2 ag 
op- 
ary H ' 
x 
rer- 
the 


686 JOURNAL OF THE AERONAUTICAL SCIENCES—DECEMBER, 


TABLE 1 
Coefficient for Symmetric Modes 
Mode & C C2 C3 
1 2.805 +0.330661 — 0.272362 —0.625031 
2 4.544 +0.585161 —0.574662 +0.396971 
3 5.459 +0.031859 —0.027602 +0.706389 
4 6.520 — 0.603591 +0.605075 +0.368340 
5 7.689 —0.574962 +0.575420 —0.411605 
TABLE 2 
Coefficient for Antisymmetric Modes 
Mode C2 C; 
1 3.105 +0.515611 — 0.473660 — 0.436791 
2 3.936 — 0.006828 +0.026366 — 0.707532 
3 5.400 —0.708691 +0.713236 — 0.064970 
4 6.908 —0.112940 +0.113940 +0.698192 
5 7.134 +0.046158 —0.045361 +0.705547 


of the sixth, so that without loss of generality one of the 
constants may be chosen arbitrarily. For convenience, 
the coefficient C, is chosen in such a manner that at 
the free end of the beam the amplitude of the bending 
modes will be equal to unity. Thus, by introducing 
another boundary condition in Eq. (5), such that, 
ats = b, 


P(s) = 1 (14) 
it can be shown that 


lcosh sin § + sinh £ cos & ) 
Y cosh § — cos & 
_ sinh sin + cosh cos — — 1 


G= 


2 cosh § — cos & (15) 
1cosh cos — sinh sin — 1 

2 cosh § — cos é 
Cs = Cr 


- In order to evaluate the above coefficients and thus 
to determine the modes of vibration from Eqs. (5), it is 
necessary to find the natural frequencies of vibration 
by solving Eq. (11) for A or é. Eg. (11) cannot be 
solved in a general form, but for values of 


a = 23° 28’ 48” | 
no = 8.8/38.8 = 0.228 (16) 
ro/b = 87/1,125 = 0.0775 


corresponding to the XB-35 wing, the first five roots 
are listed in Table 1. By substituting these values of 
& into Eqs. (15), we get the values of C;, C2, and C; for 
each of the first five symmetric modes of vibration. 
These are also entered in Table 1. 

In order to compute the modes P and Q from Fqs. 
(5), it will be more convenient to rewrite this in the 
form 


P(s)n Cin cosh £,(s/b) + C2, [sinh £,(s/b) + 
sin + cos En(s/b) (17) 
Q(s)n 2Con(En/d) cot a{cos[E,?(s/b)] + 


tan (En?) sin [En2v(s/b)] } 
where y = 1/bWmo and the subscript n refers to the 


1947 


mode numbers 1, 2, 3, 4, and 5. The above equations 
are plotted on Fig. 9. 


Antisymmetric Case 


The conditions at the origin for this case are (see 
Fig. 4) that there is no deflection, no component of 
moment reaction in the X-direction, and no rotation of 
the tangent parallel to the X-axis. The last two condi- 
tions can be recognized on Figs. 5 and 6 and can be ex- 
pressed as follows: ats = 0, P = 0, 


EI(d?P/ds?) cos a — GJ(dQ/ds) sin a = 0 (18 

(dP/ds) sin a + Qcosa = 0 ) 
By applying the boundary condition of Eqs. (6) and 
(18) to Eqs. (5) in a manner similar to that done for 


‘the symmetric case, it can be shown that the frequency 


Fic. 8. Projections of antisymmetric vibration. 
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Fic. 9. Symmetric modes of vibration. Uniform beam. 
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equation for the antisymmetric case will be 
(= Esin — sinh Ecosé Db _ 


l 
m4 2 1 
1 + cosh cos Vr 


Again in the limiting cases when a = 0 or when a = 
1/2, the above equation can be reduced to a familiar 


form. If a = 0, the frequency equation becomes cot 
Db = 0 or = 3/2z, etc.; these are the un- 
coupled torsional modes. 

For the bending modes 


cosh ky/Xb sin ky/Xb — sinh = 
or 
tanh ky/db — tan ky/dd = 0 (20) 


This equation is identical to the frequency equation 
for a ‘“‘free-hinged”’ beam given in reference 2, p. 431. 
For a = 7/2, the frequency equation becomes 


1 + cosh ky/Xb cos ky/Xb = 0 (21) 


which corresponds to the frequency equation for a beam 
with one end built in and the other end free (reference 1, 
p. 234). 

Projections of the antisymmetric modes of a vi- 
brating beam will satisfy Eqs. (20) and (21), Fig. 8, in 
the same manner as the projections of symmetric modes 
satisfy Eqs. (12) and (13), Fig. 7. 

The first five roots (é’s) of the frequency equation, 
Eq. (19), obtained using the constants of Eq. (16), are 
listed in Table 2. With the additional boundary con- 
ditions of Eq. (14)— that at the free end the amplitude 
of the bending mode equals unity—the expressions for 
the constants in Eqs. (5) for the antisymmetric case 
become: 
lcosh sin  — sinh cos & 
2 sinh £ + sin & 


_ lcosh £ cos § — sinh — sin § + 1 


1 => 


sinh + sin 
C; = 
(22) 
novo 
B, = B, cot 


Thus, for every root (&) of the frequency equation, all 
coefficients can be tabulated. The first five roots are 
listed in Table 2. 

The modes of vibration for the antisymmetric case 
can be calculated from the following equations: 


P(s)n = Ci, [cosh £,(s/d) — cos £,(s/b)] + 
Con sinh &,(s/b) + Cin sin &,(s/d) 

Q(s)n = 2C,(1/noyb) cot + 
cot cos [&y(s/b)]}} 


(23) 


The above equations are plotted on Fig. 10. 
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Fic. 10. Antisymmetric modes of vibration. Uniform beam. 


Rigid Modes 


In each of the above modes, P’s and Q’s are associated 
with definite frequencies of vibration and elastic defor- 
mations of the beam. A complete description of the 
motion of any particle must also include the motion of 
the beam as a rigid body that can be associated with 
any frequency of any mode. This type of motion 
should be determined also from Eqs. (5) if k = 1] = 
The equations 


P= Cos? + al 
Q = os + & stad 


satisfy the differential Eqs. (4) and are solutions for 
rigid modes provided the constants of integration are 
determined from the boundary conditions. 

For the symmetric case the boundary conditions are 
given by Eqs. (6) and (7), which yield 


a= 0, a= 0, o& = 0 

C3 COS a — Cssina = 0 (25) 
Substituting the above into Eqs. (24), 

P = 635 + &%, Q = c¢cota (26) 


Two independent rigid motions can be represented 
by these equations if the coefficients are properly se- 
lected. 

If c; = 0 and c = 1, then 


Py, = landQ =0 (27) 
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This is the rigid mode of a beam in which the center of 
gravity has a vertical translatory oscillation of unit 
amplitude. 

If coefficients are chosen so that « = —1 and cg = 
2/b, then 


Pw = (2s/b) — land Qoo = (28) 


This is the rigid mode of the beam corresponding to 
rotation about an axis parallel to the Y-axis passing 
through the center of gravity of the beam. This mo- 
ticn is known as the “rigid pitch.” 

For the antisymmetric case the boundary conditions 
are given by Eqs. (6) and (18), which yield 


(2/b) cot a 


qQ = 0, a = 0 
a= 0, og = 0 
sin a + = 0 
Substituting the above into Eqs. (24), 


P = —c; tana (29) 


This demonstrates that there is one antisymmetric 
rigid mode corresponding to rotation of the rigid body 
about the X-axis (roll). To show this, let c; = 1/8, 
then 


= s/band Q = —(1/b) tana (30) 


Part II.—NONUNIFORM BEAM 


Consider a nonuniform, sweptback wing as shown on 
Fig. 2. The new equations of motion will be little 
changed from those of Eqs. (1) (reference 1, p. 264) J 


+» 
>| (¥)]- (FP) 


Assume harmonic vibration, similar to that of the uni- 
form beam, then 


v 
B = exp(twt)y(s) 


where ¢(s) and ¥(s) are functions of the s-coordinate 
and determine the shape of normal modes of vibration 
of nonuniform beam. 

By partial differentiation of Eqs. (32) with respect 
to s and ¢ and substitution into Eqs. (31), the total dif- 
ferential equations of the normal modes of vibration 
may be obtained: 


d de 
wme = 0}: 


(32) 


(33) 


The solution of the above for ¢ and y will be derived 
by the Ritz method (reference 1, p. 259, and reference 
4, p. 353). 


1947 


Since the modes of vibration P and Q of the uniform 
beam satisfy the necessary boundary conditions of 
nonuniform beam, they can be used for the solution of 
g and y in the form of infinite series: 


= a,P, 


n=1 


y= anQn 


n=1 


(34) 


where the coefficients a, must be chosen in such a way 
as to satisfy the principle of minimum energy in the 
system. 

The maximum potential energy U and the maximum 
kinetic energy T in the system are: 


T = 1/2u* fo’ (me* + I,y*)ds 
which with substitution of Eqs. (34) become 


U= 1/230 me 


m=1 n=1 


(36) 
T = 1/2 mn 
where 
Aus. + 
Ban Bag’ + Bun” 
A = frEl(dP,/ds?) (d?P,,/ds?)ds\ (37) 


= ds 
ds 


The coefficients A ,,, are associated with the potential 
energy of deformation, and coefficients B,,, are -as- 
sociated with the kinetic energy of the system. The 
first are dependent on the elastic properties of the wing 
—that is, on EJ and GJ—and are independent of mass 
and its distribution; the latter depend on m and J, but 
are independent of the elastic properties. This means 
that for a given wing the integration A,,, needs to be 
performed only once for the symmetric case and only 
once for the antisymmetric case, but the integration for 
Brnn must be repeated for each different loading condi- 
tion. Integrals A,,,* and B,,,* involve the bending 
of the wing, and integrals A,,,** and B,,,** involve the 
twisting of the wing. 

In all free vibration, the amount of energy in the sys- 
tem is constant, and therefore the quantity UV — T = 0 
is always satisfied. Thus, the Eqs. (36) could be solved 
for frequency w in terms of unknown coefficients a, and 
known coefficients Am, and B,,,, but the smallest 
amount of energy involved in the vibration will occur 
when w is minimum or when the integral 


A..** = 


I= 1/2 > (Amn (38) 
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TABLE 3 
A = Matrix. Symmetric Case, Eq. (37) 
Bending Component 
394 8,305 —7,142 — 25,586 — 23,416 
8,305 15,776 1,842 — 22,042 —30,308 
At = —7,142 1,842 32,983 62,768 29,731 
— 27,586 — 22,042 62,768 163,853 119,808 
— 23,416 —30,308 29,731 119,808 137,155_] 
Twisting Component 
4 5,280 — 10,837 — 2,937 19,562 4,433 
— 10,837 41,646 — 19,450 — 50,682 84,632 
At* = — 2,937 — 19,450 37,463 —21,191 — 135,872 
19,562 — 50,682 —21,191 348,458 —28, 
| «4,433 84,632 — 135,872 — 28,806 934,104 
Composite Component 
—2,532 — 10,079 — 8,024 — 18,983 
— 2,532 57,422 — 17,608 —72,724 54,324 
A= — 10,079 — 17,608 70,446 41,577 — 106,141 
—8,024 —72,72A4 41,577 512,311 91,002 
— 18,983 54,324 — 106,141 91,002 1,071,259 


is a minimum (reference 1, p. 265). By partially dif- 
ferentiating Eq. (38) with respect to all coefficients 
a, and setting each of the resulting equations to zero, a 
series of the simultaneous equations may be obtained: 


= (Ame — Bast, = 0 (39) 
(m = 1, 2,3, ... #) 


The above series of equations have a solution for a, 
different from zero only if the determinant of these 
equations is equal to zero. This condition formulates 
the frequency equation from which the frequencies w’s 
of the various modes can be found. The solution of the 
frequency equation yields, in general, an infinite num- 
ber of values for w; Eqs. (39) have a solution only if w 
is one of the roots of the frequency equation. Therefore 
for each different w we have a different series of a,’s. 

The process of finding w’s and a,’s is laborious. Two 
general methods are available. One calls for the solu- 
tion of the frequency equation for several values of w 
and then for several solutions of series of simultaneous 
equations for a,’s. Another method, that of “matrix 
iteration,’’ is usually more convenient for this type of 
problem. The description of this method is not in the 
scope of this paper (see references 4 and 6), but the 
formulation of matrices, matrix equations, and evalua- 
tion of constants will be discussed here. 

At this point it is necessary to decide how many 
terms of the infinite series, Eqs. (29), are to be retained 
in all future work. Employment of a greater number 
of terms usually effects more accurate results. On the 
other hand, the basic assumptions of the problem do not 
justify the use of many terms, and, in addition, the 
amount of calculation necessary increases out of pro- 
portion to the number of terms retained. For these 
reasons it has been decided that the series will consist 
of the first five elastic modes of vibration of the uniform 
beam. The rigid modes, however, must also be in- 
cluded. As a result, the solution for the symmetric 


case, which has two rigid modes, must be extended over 
seven terms, while the solution for the antisymmetric 
case, which has only one rigid mode, requires six terms. 
Nevertheless, it will be shown that both cases can be 
simplified to involve the operations on square matrices 
of the fifth order. 


Symmetric Case 
The symmetric case involves the solution of the fol- 
lowing set of equations: 


(An — = 0 (40) 
(m = 1, 2, 3, 4, 5, 6, 7) 


These seven equations will contain 49 coefficients A, 49 
coefficients B, and 7 coefficients a. Both A’s and B’s 
can be arranged into square matrices of seventh order. 
Since the first two modes are rigid modes, the strain 
energy in each will be zero, and thus, from Eq. (37), 


Ai, = Ax = Am = Am = 0 (41) 
(m and n = 1, 2, 3, 4, 5, 6, 7) 


Consequently A,,, can be expressed as a square matrix 
of fifth order of elastic modes only, in which m and n = 
3, 4, 5, 6, 7. For a given airplane, the distribution of 
the elastic properties (EJ and GJ) is known. Further- 
more, d?P/ds? and dQ/ds are known at all points of the 
wing for all modes of the uniform beam; therefore, 
Amn* and Am,** can be obtained from Eq. (37) by in- 
tegration or by equivalent summation. For the air- 
plane shown in Fig. 1, calculations for A* and A** were 
performed separately, and the resultant matrices are 
shown in Table 3. It was known that the wing was 
“unusually rigid in torsion,” so it was anticipated that 
the strain energy in torsion would be small compared 
with the strain energy in bending, and therefore the 
terms A** could be neglected in comparison with the 
terms A*. No such conclusions could be justified by 
Table 3. 
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Returning to the solution of Eqs. (40) and applying 
conditions (41), the above equations can be written as: 


7 
—w? Bina, = 0 
n=1 
(42) 
—w? Bond, = 
n=2 
The remaining five equations are: 
7 
(Aus Bmn)On = 0 (43) 
n=1 
(m = 3, 4, 5, 6, 7) 


Expanding Eqs. (42) and (43), the following equations 
are obtained: 


— w > = 0 
n=3 
7 
— B =@ 


7 
n=3 
&*Bmn)d, = 0 
(m = 3, 4, 5, 6, 7) 


Coefficients a; and az can be solved explicitly from the 
first two of the above equations: 


7 
— 


(45) 
py (Bu Bin — By 
m= 
By Bu — 
(n = 3, 4, 5, 6, 7) 


so that only five unknown a,’s remain. 
Substitution of the values of a; and a, from Eqs. 
(45) into the last of Eqs. (44) results in 


7 

(Amn — = 0 (46) 
in which 


By — BuBy 
Br2BuBin "Th 
By Bn — BuBy 
(m and n = 3, 4, 5, 6, 7) 


This is also a square matrix of fifth order in which 
Brn is the matrix of the elastic modes only and can be 
evaluated by Eqs. (37). The other terms in Eq. (47) 
are functions of the kinetic energy of rigid motion of 
the body and can be evaluated as follows: 

Remembering that the subscripts in all terms satisfy 
the condition mn = nm, from Egs. (27) and (28) 


P, = 1 On = 0 
P, = 2(s/b) — 1 Q2 = (2/b) cot a 


(47) 


P,P; = 1 

P,P, = 2(s/b) — 1 
P,P, = [2(s/b) — 1]? 
Py, = P, 2:0, = 0 

Pon = [2(s/b) — 1]Pn Q:Qn = (2/6) cot aQ, 


Note: Change of notation: subscripts 0 — 1 and 
00 — 2. 


Substituting the above into Eqs. (37) the following 
coefficients are determined: 


Bu = So’ m ds 
By = fo’ m[2(s/b) — 1]ds 
Bu = So’ m[2(s/b) — 1]*ds + 
Jo’ I,(4/b?) cot? a (49) 
Bi, = So’ m P, ds 
Bon = Jo’m[2(s/b) — 1]P,ds + 
So’I,(2/b) cot a Q,ds 


where » = 3, 4, 5, 6, 7. Inasmuch as coefficients B,,, 
and, therefore, the modes of vibration are affected by 
the mass distribution along the wing, the complete sets 
of modes were calculated for three different loading 
conditions, namely: (A) minimum flying weight; 
(B) maximum landing weight; and (C) tail heavy 
flight. 

These conditions represent 46.5, 67.5, and 85.0 per 
cent of the gross weight of the airplane. Coefficients 
B* and B** were calculated by Eqs. (37) for all three 
loading conditions, and, as in the case of A* and A**, 
the magnitude of B* was of the same order as that of 
B**, showing that the kinetic energy due to twist could 
not be neglected in comparison with kinetic energy due 
to bending. 

When all entries of matrix B,,,, are collected and coef- 
ficients B,, and B2, are calculated from Eqs. (49), then 
matrix Cn, can be formulated by Eqs. (47). Having 
obtained all entries in A», and C,,,, then the matrix 
Eq. (46) can be written in the form 


(A — w*C)a = 0 (50) 


in which A and C are known square matrices of fifth 
order, w is the unknown frequency, and a is the unknown 
column matrix. 

The solution for w’s and a’s was obtained by the 
method of matrix iteration described in reference 6. A 
typical example of the results is shown in Table 4. 
Coefficients a3, a4, ds, ds, and a; and the w’s are direct 


=0 
2:0. = 0 


(4/62) cot? a } (48) 


TABLE 4 


Mode 1 2 3 4 5 
o 14.98 31.90 41.90 112.52 139.09 
a 0.57727 —3.28954 0.34112 —0.05046 0.05299 
a, 0.63921 —10.52349 0.39122 —0.76775 —0.46921 


0.22429  -—6.05682 0.13605 —0.85515 —0.92125 
0.08161 —5.31206 0.75868 0.38608 1.11555 
as 0.00769 0.16498 0.23626 —0.03380 —0.59402 
0.01076 0.49589 0.16776 —0.18406 0.35547 
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results of the iteration, while coefficients a, and a2 are 
calculated from Eqs. (45). The modes of vibration 
g and y of nonuniform beam can be obtained now by 
summation of the products of coefficients a, and modes 
of vibration P and Q of the uniform beam according to 


Eqs. (34). 


Normalized Modes 


It can be recalled that, in the calculation of the modes 
of vibration of the uniform beam, it was necessary to 
introduce an arbitrary constant that was chosen in such 
a way as to give unit amplitude at the free end. A 
similar requirement of an arbitrary constant is neces- 
sary in the calculation of the modes of the nonuniform 
beam. This constant appears in each iteration as 
unity for the coefficient a3 (see Table 4). This means 
that all other a,,’s and, consequently, the amplitudes of 
the modes (see Eqs. 34) are measured in terms of units 
of a3. Since this is an arbitrary scale, it can be changed 
to some other more convenient standard. It will be 
illustrative to compare the various modes of vibration 
on the basis of the amount of energy in each. Since the 
kinetic energy is as given by Eq. (35), then, with the 
stipulation that the amount of energy in each new mode 
¢, and y, is equal to unity, the above conditions can be 
expressed as 


fo” (me? + = 1 


where the subscript j refers to the mode number. 
It is easy to deduce by substitution that ¢,; and Y, are 
related to ¢; and y; by the expressions 


= ¢;/N,and J; = ¥,/N; 


(51) 


(52) 
where 
N? = 1/207 fo” (me? + (53) 


is a constant for each mode and is known as the “‘nor- 
malizing factor.” 

Symmetric normalized modes of vibration of a non- 
uniform beam for three different loading conditions are 
shown on Fig. 11. 


Antisymmetric Case 


In the antisymmetric case there is one mode of rigid 
motion, and, as a result, for five elastic modes it is nec- 
essary to solve six simultaneous equations similar in 
form to Eqs. (40). 


6 
n=1i1 


(m = 1, 2, 3, 4, 5, 6) 


(54) 


The analysis of this case is parallel to that of sym- 
metric modes; the discussions are similar. 

For the first rigid mode it can be shown by means of 
Eq. (37) that 


(55) 


Ai, Am 
(m and n = 1, 


0 
2, 3, 4, 5, 6) 


BENDING COMPONENT Q(s) TWISTING COMPONENT 


CENTER LINE Tip CENTER LINE TIP 
OF BEAM OF BEAM OF BEAM OF BEAM 
LOADING 
----M 
—---D 


MODE ~ 


Nonuniform beam. 


Fic. 11. Symmetric modes of vibration. 


Therefore A,,, can be expressed as a square matrix of 
fifth order of elastic modes. With the above condition 
the first equation of Eqs. (54) reduces to 


6 
—w? = 0 (56) 
while the remaining five are . 
6 
(A mn Bin )On =0 (57) 
n=1 


(m = 2, 3, 4, 5, 6) 


Expanding Eqs. (56) and (57), Eqs. (58) may be ob- 
tained 


6 
n=2 


(58) 
+ (A mn = 0 
n=2 
(m and n = 2, 3, 4, 5, 6) 
from which 
6 
a=- Bintn/ Bu (59) 


Consequently, only five unknown a,’s remain. 
Substitution of the value of a; from Eq. (59) into the 
second of Eqs. (58) results in the equation 


| 
MODE — 
“A / 
/ 
ed by MODE ~ 
te sets a _ 
vading | | / | 
eight; 4T MODE ~ 
heavy 
cients 
three 
| At, 
hat of 
could 
due | 
| coef- 
, then 
aving 
natrix 
(50) 
fifth 
y the 
6. A | 
4. 
direct 
; a 
).09 
95299 
46921 
92125 
5940 a 
35547 | 


692 JOURNAL OF THE AERONAUTICAL SCIENCES—DECEMBER, 1947 


BENDING COMPONENT @(s) TWISTING COMPONENT \y(s) 


CENTER LINE TIP CENTER LINE TIP 
OF BEAM OF BEAM OF BEAM OF BEAM 


LOADING 


4, 


RO MODE 


MODE 


‘usm mopE—/ 


Fic. 12. Antisymmetric modes of vibration. Nonuniform beam. 


6 
(Ams mn) On = 0 (60) 
n= 
(m and n = 2, 3, 4, 5, 6) 
This is also a square matrix of fifth order in which B,,, 
is the matrix of the elastic modes evaluated from Eqs. 


(37). The other terms in Eq. (61) are functions of the 
kinetic energy of motion, as a rigid body, and can be 


found as follows: From Eq. (30) 
= s/b, = —(1/d) tan a 


or 


P,P; = (s/b)? 0:02 = (1/b?) tan? \ (61) 
PiP, = (s/b)Pn  QiQ, = —(1/6) tan a Q, 


Substituting the above into Eq. (37) the following coef- 
ficients are determined 


So’m(s/b)*ds + So'l,(1/b?) tan? a ds (08 
— tan aQ,dsf ©) 


Thus, Eq. (54) can be stated in exactly the same form 
of matrix equation as Eq. (50) for the symmetric case, 
The solution of this equation would finally lead to the 
normalized modes of vibration for antisymmetric case 
shown in Fig. 12. 

In conclusion, it may be noted that the above analysis 
in its original form involved approximately one quarter 
of a million multiplications of two numbers of six sig- 
nificant figures. This task was performed with only a 
small expenditure of effort and time by means of the 
punched card system of the International Business 
Machines. 
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Bin 
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Heat-Resistant Alloys for Use in 
Jet Propulsion Engines 


J. W. FREEMAN,* E. E. REYNOLDS, ano A. E. WHITEt 
Unwersity of Michigan 


INTRODUCTION 


pense OPERATING at high temperatures are ex- 

pected to deform progressively under the stresses 
encountered and to deteriorate progressively from 
corrosion and structural instability. Their suitability 
for any given application depends on the rate of deforma- 
tion or deterioration being sufficiently slow so that the 
service loads can be carried for the expected life of the 
unit without failure. Other requirements are normal 
in that the material must be commercially available in 
the necessary form, can be fabricated, and will not be 
excessive in cost. 

Aircraft propulsion engines are, as yet, a relatively 
short-life application. Operating stresses and tempera- 
tures are kept high for efficiency and to save weight. 
Corrosion is not too difficult a problem in most cases 
because the materials used, on the basis of strength con- 
siderations, generally have sufficient corrosion resist- 
ance. 

Basically, the aircraft engine is a short-life high- 
stress application. Uncertainties due to extrapolation 
of data, estimation of structural stability, and corrosion 
effects are relatively minor because the service life is as 
yet no more than the time periods customarily used in 
laboratory evaluation of high-temperature materials. 
The main problem has been to develop alloys with high 
strength at high temperatures for relatively short time 
periods and to be able to produce them in quantity with 
high-quality and uniform properties. 

The fabrication requirements, as well as other con- 
siderations, such as heat-shock resistance, have limited 
the materials used to metallic alloys. Traditionally, 
the development of the gas turbine has been retarded 
by the low strength of available high-temperature 
alloys. Up until 1941, presumably these considerations 
had limited the gas-turbine development in this country 
to an experimental aircraft turbosupercharger, a Diesel 
turbosupercharger, and low-temperature gas turbines 
operating on waste process gases. 

The urgent requirements for aircraft turbosuper- 


chargers and jet engines during the last war initiated 


vast research programs relating to the metallurgy of 


heat-resisting alloys. Many new alloys were developed 


Presented at the National Aircraft Propulsion Meeting, 
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with strengths much higher than had previously been 
available. In addition, valuable design experience was 
obtained which makes possible more efficient utiliza- 
tion of the available materials and greater confidence in 
the basis for design. 

In 1941, the background for the development of al- 
loys with high strength at elevated temperatures was 
based almost entirely on experience in the steam-power, 
oil, chemical, and furnace construction industries. 
These were almost invariably long-life applications at 
relatively low stresses, with practically no dynamic 
loading. The alloys used included such materials as 
carbon steel, !/, per cent Mo steels with Cr contents up 
to 9 per cent, and stainless steels of the 12Cr, 18Cr, 
27Cr, 18Cr-8Ni, 25Cr-12Ni, 25Cr-20Ni, and 35Ni- 
15Cr types. Some experience with 65Ni-15Cr, 80Ni-20 
Cr, and Inconel alloys was available. 

Even at that time some developmental work had 
been done which had led to the availability of stronger 
alloys at high temperatures, such as the Universal- 
Cyclops’ alloys 17W and 19-9WMo used in the early 
turbosuperchargers and in the waste gas turbines. The 
Timken Roller Bearing Company had been experi- 
menting with the then phenomenally high-strength 
16Cr-25Ni-6Mo alloy. Hastalloy B was known but had 
been advocated primarily as a special corrosion resist- 
ance alloy. K42B was being developed at Westing- 
house. : 

The limited experience with heat-resisting materials 
in gas turbines was contradictory to most experience, 
in that the aircraft turbosuperchargers had a short- 
life requirement with high stresses. On the other hand, 
the Diesel turbosupercharger and waste-gas-turbine 
experience was not particularly applicable because these 
units were designed for long life at relatively low stresses 
and temperatures by the same considerations used for 
design in steam-power and oil-refinery practice. 

Some of the metallurgists working with the alloys 
were beginning to realize at that time that precipitation 
reactions of complex carbides and intermetallic com- 
pounds were probably a controlling factor in the high- 
temperature strength of the more complex alloys. This 
was little better than a hazy realization of the probable 
explanation. Other general metallurgical considera- 
tions were available for use: refractory metals with 
stable structures were needed; austenitic base alloys 
were inherently stronger at high temperatures than the 
ferritic alloys; molybdenum, tungsten, columbium, and 
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titanium would improve the strength of steel at elevated 
temperatures. 

With these considerations in mind, a number of 
metallurgical organizations started out to develop alloys 
that wouid have better strength at higher temperatures 
than known alloys. A large number of compositions 
were made up on an empirical basis and tested. A few 
of these exhibited unusually high strengths in the tem- 
perature range from 1,200° to 1,500°F.! Those that 
exhibited such unusual properties received more atten- 
tion and became known as “Super Alloys.”” Those 
well known today include Timken Alloy, Hastalloy B, 
Vitallium, 19-9DL, Low-Carbon N155, S495, S590, 
$816, Inconel X, 422-19, K-42B, Refractaloy 26, 
etc. This list is by no means complete but men- 
tions some of those that have been used to some 
extent. 

The struggle to develop large-scale production of 
these alloys has been a tedious and expensive affair. In 
the wrought alloys, troubles were encountered in ob- 
. taining good quality; they were difficult to hot-work 
by virtue of their high strength at elevated tempera- 
tures; machining was difficult; inspection was difficult 
and critical; their strength and ductility at high tem- 
peratures were erratic; and it was difficult to uniformly 
reproduce commercially the high strength of the experi- 
mental alloys. In some instances increases in strength 
were offset by increases in density. 

The materials actually used for each engine part were 
selected on the basis of the best compromise of all the 
metallurgical requirements, because any one material 
usually does not best fulfill all needs. Each engine part 
operating at high temperatures involves special con- 
siderations unique to its particular operating conditions, 
function, and design. 


Rotor Discs 


Research metallurgists have been asked to provide 
materials that will meet the service requirements of 
critical conditions at two locations in the discs: the 
high stresses at the hub due to rotation of the rotor; 
and the high stresses at the rim due to the pull of the 
blades on their fastenings. In the turbines usually used 
as part of aircraft engines, the hub is cool, but the de- 
sign is such that the stresses are high. The operating 
temperature of the rim has usually been estimated 
to be about 1,200°F. 

Yield strength at room temperature has generally 
been the criterion given to the metallurgist for hub 
service. There is no exact value applicable to all de- 
signs. In general, however, yield strengths as high as 
80,000 to 100,000 Ibs. per sq.in. for 0.02 per cent offset 
have been required in this country. This requirement 
has been high not only for service stresses but perhaps 
primarily to withstand a high overspeed test for qual- 
ity. Good ductility is needed to allow adjustment to 
stress concentrations. 
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Creep strength is the major consideration for rim 
service. Exact operating stresses are not generally 
available. Rupture strength has been widely used as a 
criterion of service performance with the tendency to 
force the stress for rupture in 100 and 1,000 hours as 
high as possible. In order to receive consideration, a 
new alloy has had to have rupture strength at 1,200° 
F. in excess of 40,000 Ibs. per sq.in. for fracture in 
1,000 hours. Elongation in the rupture test is also 
considered important with the tendency to require 
that it be kept higher than 5 per cent to prevent pre- 
mature failures at stress concentrations. 

Such properties aS corrosion resistance, fatigue 
strength, damping capacity, and heat-shock resist- 
ance have not been major considerations in disc 
material development. 

So far it has been possible to meet the requirements 
for discs only by the use of forgings. The ailoys used 
must also be machinable. Volume production makes 
it extremely desirable to be able to weld buckets and 
shafts to the discs. 

Metallurgically, the operating temperature of 1,200° 
F. or higher, together with the required yield and rup- 
ture strengths, has imposed severe restrictions. Heat- 
treatable ferritic alloys cannot meet the rupture 
strength requirement. Standard austenitic steels also 
have too low rupture strengths. The high-creep- 
strength austenitic alloys, which contain molybdenum, 
tungsten, columbium, titanium, cobalt, etc., must be 
cold-worked in order to develop suitable low-temperature 
yield strengths. The theoretically ideal age-harden- 
able austenitic base alloys have not been adapted to the 
service, one of the reasons being their low elongation and 
reduction of area in the rupture test. One solution 
to the problem has been to weld an austenitic steel 
rim on a hardened ferritic steel center thus forming a 
composite disc. 

The highly alloyed austenitic base Timken Alloy and 
19-9DL have been most widely used in this country 
and G18B in England. Their properties at high tem- 
peratures have certainly been variable, and yet the 
service records of discs have been good. The metal- 
lurgist is therefore inclined to regard the customary 
quoted conditions of disc service somewhat critically. 


There seems to be reason to believe, however, that as 
experience is acquired it will be possible to design tur- 
bines in which the upper temperature limits of the 
rotor discs will be considerably restricted to take ad- 
vantage of the lower cost, high quality, and good fab- 
rication characteristics of heat-treatable ferritic alloys 
as has been done to some extent in England and Ger- 
many. The certainty of increased operating tem- 
peratures, particularly in turbines for driving pro- 
pellers, indicates a strong need for both the austenitic 
and age-hardenable alloys. For these reasons basic 
metallurgical principles of these three types of alloys 
will be reviewed. 
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HEAT-RESISTANT ALLOYS FOR USE IN JET ENGINES 695 
an TABLE 1 
Dise Alloys 
Chemical Composition, Per Cent ~ 
Alloy | & Si Mn Cr Ni Co Mo WwW Cb Ti N Cu 
Austenitic Base Alloys 
Timken 0.10 0.75 1.00 16 25 6 0.15 
19-9DL 0.30 0.60 1.00 19 9 m. 1°35 1 0.4 0.4 
Low-Carbon N155 0.15 0.40 1.50 20 20 20 3 2 1 da 0.15 
$590 0.45 0.50 1.00 20 20 20 4 4 4 a 
S816 0.45 0.50 0.75 20 20 45 4 4 4 < r 
EME 0.15 0.50 0.40 19 11 ry 3.5 1.3 0.10 ‘ 
CSA 0.40 0.30 4.00 18 5 ae 1.5 1.5 0.6 te ; 
*G18B 0.40 1.00 0.80 13 13 10 2.0 2.5 3.0 bf 4 
*Rex 78 0.08 14 18 4 0.6 4 
Age-Hardenable Alloys 
Fe Al 
Discalloy 0.05 0.70 0.70 13 25 Bd 3 55 - 1.8 0.20 
Inconel X 0.05 0.40 0.50 15 73 <s ee 7.0 1.0 2.5 0.70 
* Alloys developed in England. 
TABLE 2 


Selected Properties of Austenitic Base Alloys! 


Room Temp. 
0.02 Per Cent Offset Stress for Rupture in 1,000 Hours 


Yield Strength Lbs. per Sq.In- 
Alloy Treatment (Lbs. per Sq.In.) 1,200°F. 1,350°F. 1,500°F. 
Timken Solution-treated 47,000 35,000 ae 8,000 
Hot-cold worked 119,000 42,000 19,000 et 
19-9DL Solution-treated 26,500 34,500 sagt 10,000 
Hot-cold worked 100,000 45,000 20,000 ead 
Low-Carbon N155 Solution-treated 42,000 37,000 22,500 13,000 
Hot-cold worked 135,000 48,000 27,000 reegepe 
$590 Solution-treated + aged 50,000 42,000 25,000 15,000 
$816 Solution-treated + aged 68,000 45,000 30,000 18,000 
EME Hot-cold worked 100,000 40,000 
CSA Solution-treated 50,000 37,000 12,000 ‘wwe 
Hot-cold worked 90,000 32,000 
*G18B* Solution-treated 34,000 t 37,000 19,000 10,000 
*Rex 78* Solution-treated + aged 48,000 Not available 
* Alloys developed in England. 
+ 0.2 per cent proof stress. 


Hardenable Ferritic Steels 


The ferritic steels have not received much considera- 
tion in this country. Consequently, there are rela- 
tively little data of the type usually supplied for jet- 
engine materials selection. There is, however, no 
question but that such alloys would be the most suit- 
able if operating temperatures could be limited to 
about 1,000°F. Their strength falls off rapidly because 
of tempering with further increase in temperature. 
There would be no difficulty in developing suitable low- 
temperature physical properties. 

Metallurgically, the best steels would probably be of 
the Cr-Mo-Va type with 0.2-0.4 per cent carbon. 
Experience indicates that the best high-temperature 
properties would be obtained by normalizing (air 
cooling) from 1,700-1,750°F. and tempering, rather 
than oil quenching and tempering back to the desired 
physical properties. 

Such steels would be relatively easy to produce uni- 
formly with good quality and could be available in 
large quantities, low in cost, and relatively easy to 
fabricate into turbine discs. Their low expansion 


would decrease thermal stresses. 
are so pronounced that engine designers will certainly 
give careful consideration to the possibilities of cooling 
and other means to permit their utilization in engines 
without exceeding their limiting service temperature 
of about 1,000°F. 


Approximate Properties of Possible Ferritic Disc Materials 
Temperature, °F. 75 800 900 §61,000 1,100 
0.2 per cent offset yield strength 110,000 68,000 65,000 62,000 


Stress for rupture in 1,000 hours 70,000 60,000 50,000 32,000 
1 per cent creep in 1,000 hours ys 55,000 45,000 25,000 
0.1 per cent creep in 1,000 hours 48,000 30,000 15,000 

The importance of heat-treatment, however, is 
shown by the fact that this type of steel, when in the 
oil-quenched and drawn condition under a stress of 
50,000 Ibs. per sq.in. at 1,000°F., fractured in 250 
hours; whereas, when in a normalized and drawn con- 


dition, it did not fracture in 1,000 hours. 


Austenitic Base Alloys 


This class of alloys includes those materials that 
contain sufficient nickel and chromium to develop 
stable austenite; have various combinations of molyb- 


These advantages 
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denum, tungsten, cobalt, columbium, and titanium 
added to enhance the strength at high temperature; 
and do not harden substantially without cold-work. 
Typical examples of alloys of this type are shown in 
Tables 1 and 2. 

The principles involved in the metallurgy of these 
alloys are not well understood. Such generalities as 
are given below are limited because there are many 
exceptions for specific alloys, and in many cases they 
are based on limited experience. 

(1) Yield strengths at room temperature are limited 
to values of 30,000 to possibly 60,000 Ibs. per sq.in. 
unless hot-cold work is applied. Reductions of 15 to 
20 per cent at temperatures between 1,500°and 1,200°F. 
will raise the yield strengths to above 100,000 Ibs 
per sq.in. Even after this hot-cold work, ductility 
characteristics are usually adequate. 

(2) The hot-cold work increases the strength at 
temperatures up to 1,350°F. under high stresses and 
for the time periods usually considered in aircraft tur- 
bines. 

(3) The combination of properties obtained in any 
one material at any one temperature depends on the 
prior thermal treatment, the temperature of hot-cold 
work, and the degree of reduction. 

(4) A solution treatment prior to hot-cold work 
usually increases strength at high temperatures but 
reduces the elongation and reduction of area in the 
rupture test. This effect is greater, the higher the 
solution-treating temperature. 

(5) Comparisons between the as-forged and the 
solution-treated condition are apt to be misleading. 
A high finishing temperature during hot-working may 
be equivalent to an effective solution treatment. Sub- 
sequent tests may give altogether different resylts from 
those obtained from a material with a lower finishing 
temperature which promotes prior precipitation and 
hot-cold work. 

(6) The best combination of strength and duc- 
tility is occasionally found in as-forged materials with 
low finishing temperatures.? This process has appar- 
ently been difficult or impossible to control and dupli- 
cate; consequently, production of discs has been limited 
to reheating to temperatures between 1,200° and 1,500° 
F. for hot-cold working. 

(7) Fabrication procedures that involve repeated 
reheating to temperatures between 1,700° and about 
1,900°F. after hot-cold work may reduce the strength 
to a pronounced extent. This may lower the strength 
below the inherent strength of the unalloyed aus- 
tenite. 

(8) For service at temperatures above 1,350°F., for 
other than extremely short time periods, the best prop- 
erties are obtained by a combination of solution treat- 
ment and aging. The considerations involved in this 
type of service overlap into the field of present-day 
turbine blades, and they will be considered together in 
the next section. 
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(9) The basic metallurgical reasons for the various 
effects are not definitely known. A working theory 
has been developed on the basis that the strength at 
either low or high temperatures is dependent upon 
strain hardening and the composition, size, and dis- 
tribution of carbides, nitrides, and intermetallic com- 
pounds. The relative importance of these effects is not 
known except that the changes in high-temperature 
properties of plain chromium-nickel austenitic alloys, 
as a result of hot-cold work, are negligible in comparison 
to changes in alloys containing Mo, W, Cb, and Ti, 
which are strong formers of carbide and other com- 
pounds. 

Presumably, the drop in properties at temperatures 
above about 1,350°F. as a result of hot-cold work is 
due to structural instability. Relief of strain-hard- 
ening and agglomeration of precipitated particles to a 
low-strength dispersion are the principal effects. In 
some cases, transformation to a weaker structure may 
be a contributing cause. 

(10) The wide range in strengths for any one alloy, 
depending on processing procedure and _heat-treat- 
ment, has made it difficult to determine the specific 
effects of alloying elements and to arrive at optimum 
compositions. 

Alloying elements like Cr, Mo, W, Ti, and Cb tend 
to promote a ferritic structure. The proportions of 
these elements must be so balanced with austenite 
formers like Ni, C, Co, Mn, and N that the structure 
remains austenitic, or very nearly so. There are 
two primary reasons for requiring a stable austenitic 
structure. One is that ferrite is weak at high tem- 
peratures and, in two-phase structures, during service 
at high temperatures, tends to become brittle by carbide 
precipitation and by transformation to sigma phase. 
The other is that a mixed structure may be difficult to 
hot-work. 

There is considerable evidence to indicate that a 
borderline composition, such as 19-9DL and CSA 
alloys, with the presence of small amounts of ferrite, 
may have outstanding properties. There is no way 
of proving, however, whether this is the basic cause or 
whether such analyses contain critical amounts of 
strengthening elements like Mo, W, Cb, or Ti. 

Chromium, nickel, and cobalt can vary over wide 
limits with relatively little effect on the strength at 
high temperatures. The elements molybdenum and 
tungsten, either alone or together, are essential to the 
strength at high temperatures. Columbium or titanium 
appear in most of the better-known alloys, not as carbide 
stabilizers for the prevention of sensitivity to corrosion 
but as strengthening agents. In both cases the 
formation of carbides and intermetallic compounds 
seems to be the most important reason for their strength- 
ening effect. Carbon and nitrogen do not have readily 
definable roles, other than their austenite stabilizing 
effect. Optimum amounts of carbon seem to be neces- 
sary for good properties for any one combination of 
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HEAT-RESISTANT 


ALLOYS 


elements. Nitrogen also has beneficial effects, an im- 
portant one of which is its improvement in quality. 
The high nitrogen contents are required in such alloys 
as Timken Alloy and Low-Carbon N155. Even less 
is known about manganese and silicon. Apparently, 
however, the primary role of these elements is to pro- 
mote a higher quality product. These amounts have 
therefore been fixed at the levels that steel-making 
practice has shown to be best for high yields of sound 
material. 

(11) The reasons for the variations in properties 
of any one material with thermal treatment and proc- 
essing procedure are complex. Strength at high 
temperatures is certainly a function of the size and 
dispersion of precipitated phases. The heat-treatment 
and processing conditions are believed to control 
properties through their effects on the precipitation 
reactions. The complexity of the relationships is 
believed to be the reason why the relative effects of 
chemical composition in alloys has defied systematic 
correlation. 

In addition to these strength considerations the 
alloys have certain other characteristics that are im- 
portant to their use as rotor discs: 


(1) The cost is high. 

(2) They are difficult to hot-work and to hot-cold 
work by virtue of their high strength at elevated tem- 
peratures. 

(3) They are usually hot-cold worked rather than 
cold-worked at room temperature because the reduc- 
tion in strength at 1,200° to 1,500°F. is an important 
feature in reducing the difficulty of cold-working the 
material. An even more important reason is that the 
properties are considerably enhanced by working at 
these temperatures where precipitation occurs. 

(4) Ingot segregations and the flow of the metal 
inherent to the usual upsetting operation in forming 
discs from billets introduce quality problems. The 
disc centers in particular are prone to have low ductility 
and strength, especially in large sizes. 

(5) Expensive surface preparations and extreme 
care are required during hot-working to prevent the 
formation of cracks and other flaws. 

(6) They are difficult to machine and weld. Weld- 
ing is particularly important to any proposed high- 
volume production. 


Age-Hardenable Alloys 


The age-hardenable alloys offered for gas-turbine 
service provide both high yield strength at room tem- 
perature and high rupture strength at elevated tem- 
peratures. In addition, they may be machined while 
in the soft solution-treated condition prior to aging. 

Typical examples of alloys of this type are included in 
Tables 1 and 3. Other age-hardenable alloys de- 
scribed later also have attractive properties but have 
not been advocated for discs. 
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TABLE 3 
Selected Properties of Age-Hardenable Disc Alloys 


Stress for Rupture 
Room Temp. in 1,000 Hours 


0.2 Per Cent ———Lbs. per Sq.In.—— 


Offset Yield 1,200° 1,350° 1,500° 
Alloy Treatment Strength F. F. F. 
Discalloy* Solution-treated 
ag 100,000 38,000 eee 
Inconel X’ Solution-treated 
+ aged 102,000 69,000 42,000 18,000 


The heat-resistant age-hardenable alloys have an 
austenitic base and usually contain titanium and alumi- 
num to furnish the precipitating phase. Columbium 
and molybdenum also act as aging elements in some 
alloys. Nickel and chromium are basic elements in 
the matrix and may contain various amounts of iron 
or be essentially iron free. Cobalt may be added with 
important beneficial effects. 

A comprehensive discussion of the metallurgical 
principles involved in these alloys was given by Howard 
Scott and R. B. Gordon in a paper before the A.S.M.E. 
last December.* They stated that various factors 
influencing strength at high temperatures were found 
to be as important as composition in determining prop- 
erties. Not until these factors were brought under 
close control could the effects of composition variation 
be determined. Among the factors that influence the 
properties are prior treatment, solution-treating tem- 


- perature and time, cooling rate from solution treating, 


and aging temperature and time. The specific effects 
of all of these depend upon the service temperature, 
life, and permissible deformation. The magnitude of 
the effect on properties, of course, varies to a con- 
siderable extent. With all of these variables normal to 
any one combination of alloying elements, it is not sur- 
prising that there has been some difficulty in estab- 
lishing preferred analyses and ensuring the desired 
properties in rotor discs under commercial production 
conditions. 

The complex nature of the relationships of properties 
to treatment, composition, and service conditions makes 
it impossible to generalize or to discuss these relation- 
ships fully in a paper of this type. There are, however, 
a few guiding principles that may be useful to engine 
designers. ‘ 

(1) The combinations of properties are attractive 
and should receive careful consideration. The suc- 
cessful application, however, will depend on the careful 
control of fabrication and heat-treatment, with special 
attention to establishing the correct treatments for the 
conditions of each application. 

(2) Each alloy has certain peculiarities that in- 
fluence its suitability for various service conditions. 
Overaging above certain temperatures with rapid loss 
in strength is a most important such characteristic. 
The titanium-aluminum hardened alloys are sufficiently 
stable for service at temperatures ranging for different 
alloys up to about 1,400°F. 
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TABLE 4 
Blade Alloys 
Chemical Composition, Per Cent 
Alloy _S Si Mn Cr Ni Co Mo WwW Cb | Al 
Austenitic Base Alloy 

N 
Low-Carbon N155 0.15 0.4 1.5 20 20 20 3 2 1 0.15 
$590 0.45 0.5 1.0 20 20 20 4 4 4 
$816 : 0.45 0.5 0.7 20 20 45 4 4 4 
Hastalloy B 0.05 0.7 0.7 65 

Age-Hardenable Alloys 

Fe 
K42B 0.05 0.7 0.7 18 42 22 rv 2 di 2.2 0.2 14 
Refractalloy 26 0.05 0.7 0.7 18 37 20 3 a a 2.8 0.2 18 
Inconel X 0.05 0.4 0.5 15 73 eu oe aye 1.0 2.5 0.7 7 
Nimonic 80* 0.04 0.5 0.5 21 74 2.4 0.6 43 
Refractalloy 70 0.05 0.2 2.0 20 20 30 8 4 ae ge “ee 15 

Cu 
Rex 78 f 0.08 14 18 4 0.6 4 

Precision Cast Alloys 

Fe 
Vitallium 0.20 0.6 1.0 25 3 61 5.5 2 
61 0.45 0.6 0.6 25 70 5 1 
422-19 0.55 0.6 0.3 25 10 55 6 1 
6059 0.40 0.8 0.2 25 35 35 5 1 
X-40 0.50 0.7 0.6 25 10 55 7 0.5 
$816 0.45 0.5 0.7 20 20 45 4 4 4 3 


* Analysis reported in U.S.A. work. 
English alloy. 


(3) Most of the precipitation-hardened alloys have 
low elongation and reduction of area in rupture tests. 
This may or may not be a serious limiting considera- 
tion, depending on the design and service conditions. 

(4) Low elongation and reduction of area in the 
rupture test is not necessarily an indication of low 
ductility under sudden overloads during service. In 
fact, most of the alloys will bend and stretch, rather 
than fracture, when sudden impact loads are applied 
during service. 


Rotor BLADES 


The main criteria of blade service are generally be- 
lieved to be creep strength and hot ductility. There 
is considerable indication that fatigue strength at 
elevated temperatures is likewise important. Such 
effects as thermal stresses due to nonuniform heating, 
heat shock, corrosion, erosion, and impact resistance 
have received less consideration in alloy development. 
In service, the blade is believed to be subjected to 
tensile stresses because of rotation, steady gas bending 
stresses, and alternating gas bending stresses as the 
blades pass nozzles. Excessive vibrations at natural 
frequencies of the blades may be a problem. Service 
temperatures range upward from about 1,300°F., with 
the upper limit desired for efficiency being very high. 

The stresses and temperatures are so high in blades 
that only the austenitic-type alloys have sufficient 
strength. The possibility of cooling blades to permit 
the use of low-alloy material exists, but this method 
has not been widely used in this country, and this 
discussion will be confined to alloys used in blades with- 
out cooling. Three types of alloys are used: (1) aus- 


tenitic base alloys that do not harden extensively with- 
out cold-work; (2) precipitation-hardening austenitic 
alloys; and (3) precision-cast cobalt-base alloys. 
Typical alloys that are being used or have received 
serious consideration on a developmental basis are 
listed in Tables 4 and 5. These tables do not propose 
to cover all possible alloys but are intended to be il- 
lustrative. 

Before taking up the metallurgical characteristics of 
each type of alloy, certain characteristics of blade 
service should be reviewed. Creep characteristics are 


TABLE 5 
Selected Rupture Strengths for Blade Alloys 


Stress for Rupture in 1,000 Hours 


. per Sq.In. 
Alloy 1,350° F. 1,500°F. 1,600°F. 
*Austenitic Base Alloys! 
Low-Carbon N155 23,000 12,500 9,000 
$590 25,000 15,000 9,000 
$816 30,000 20,000 9,500 
Hastalloy B 10,000 
*Age-Hardenable Alloys 
K42B 27,000 11,000 
Refractalloy 26 18,000 
Inconel X? 42,000 18,000 
Nimonic 80! 10,000 
Refractalloy 70° 24,000 15,000 10,000 
Rex 78 Not available 
?Precision-Cast Alloys! 
Vitallium 25,000 15,000 13,000 
1 36,000 ,000 12,000 
422-19 ,000 22,000 15,000 
X-40 33,000 24,000 18,000 
$816 29,000 21,000 15,000 
6059 30,500 19,000 12,000 


_* Data for wrought alloys on solution-treated and aged mate- 


+ Cast materials were cast and aged. 
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of primary importance in preventing failure by ex- 
tension to cause rubbing on the casing. The role of 
fatigue strength is not clear. The relative importance 
of design to avoid undue vibratory stresses and the 
inherent fatigue strength of the alloys have not been 
completely established. Also, there is considerable 
doubt about the proper way to run fatigue tests to 
measure the performance of an alloy in blade service. 

The importance of ductility in the rupture tests is an 
unsettled issue. The restriction to a minimum of about 
5 per cent elongation imposed by many engineers in 
this country may be an unnecessary handicap to many 
otherwise good alloys. The apparently excellent ex- 
perience of the English with Nimonic 80, in spite of 
its low elongation in the rupture test, throws con- 
siderable doubt on the high rupture-test ductility re- 
quirement. 

The production of blades presents difficulties in- 
herent in the types of alloys that have suitable prop- 
erties. They are difficult to hot-work. None of the 
alloys is readily machinable. The turbine blade is an 
ideal product for the precision casting process; yet, as 
blade sizes increase, it becomes increasingly difficult to 
make sound castings with good structure. 

In wrought blades the properties at high temperatures 
are remarkably dependent on thermal history, hot- 
working conditions, and heat-treatment. The pouring 
temperature, mold temperature, blade size and shape, 
and heat-treatment all may have a pronounced effect 
on the properties of cast blades. The relative effects 
of the different variables also depend on the service 
conditions. 


Austenitic Alloys 


All of the principles covered in the discussion of this 
type of alloy in the disc section apply to the blade metal- 
lurgy. There, are, however, certain concepts unique 
to blades. 

(1) One major difference is that high yield strengths 
are not necessary in blades. The need to use or avoid 
hot-cold work then depends upon its effect on high 
temperature strength. In general, hot-cold work ap- 
pears to be detrimental to high temperature strength 
at temperatures above 1,350°F., except for short time 
periods. Because most blades operate in this tempera- 
ture range, it has not been used. 

(2) In general, the best properties, as measured by 
rupture and creep tests, are developed by solution and 
aging treatments. As-forged blades may have variable 
properties. They will be adversely affected by hot- 
cold work induced by low forging temperatures and by 
agglomeration of precipitated particles. 

(3) The higher the solution-treating temperature, 
the higher the creep and rupture strengths at tempera- 
tures above about 1,350°F. In general, rupture test 


ductility decreases with solution-treating tempera- 
ture. 
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(4) Aging at about 1,350°F. prior to testing gen- 
erally gives the highest rupture and creep strengths 
over the temperature range from 1,200° to 1,500°F. 

(5) As further information is obtained, however, it 
will probably be found that the aging temperature 
should be adjusted to the service temperature. 

(6) An extremely high solution-treating tempera- 
ture may develop extremely high rupture and creep 
strengths. Other mechanical properties, such as fa- 
tigue strength, impact resistance, ductility, and stress- 
concentration brittleness, may offset this advantage. 
Generally, such adverse changes are accompanied by 
grain-coarsening. 

(7) Coarse grain size in itself may be the cause of, 
or may simply be an indication of other changes that 
cause, high strength and poor mechanical properties. 
There is evidence to support both views. It is prob- 
able, however, that grain size in itself is not important to 
strength at high temperatures, although it is important 
to mechanical properties. In fact, it may be that the 
high columbium alloys, such as S590 and S816, may be 
particularly good because their resistance to grain 
growth permits high solution-treating temperatures 
without excessive grain growth and the development of 
adverse mechanical properties. 

Hastalloy B and S816 have been the most widely 
used of this type of alloy. Many considerations have 
entered into their choice. Full recognition of the 
variables present would suggest that either they have 
inherently good properties at high temperatures, which 
give adequate minimum strength in spite of processing 
variables, or they respond well to established processing 
procedures. It may be that the rupture test is not too 
good an index of actual service requirement. One 
important consideration has been that their forging 
characteristics have been better than those of other 
alloys considered. 


Age-Hardenable Alloys 


The general metallurgical considerations covered 
under disc materials also apply to blading application 
and will not be covered further except for some con- 
ditions peculiar to blades. 

(1) The ease of obtaining high yield strengths in 
this type of alloy is not nearly so potentially important 
as in disc service, because high yield strength is not 
required for blades. 

(2) Low ductility in the rupture test may not be 
so important in blades as in discs. 

(3) Most of these alloys have high ductility when 
suddenly overloaded in service, even though the rup- 
ture test ductility is low. 

(4) Service temperatures for blades may be at or 
above the temperatures where titanium-aluminum 
age-hardenable alloys overage and rapidly lose strength. 

K42B and Refractalloy 26 have been used in this 
country. Nimonic 80 has been successful in England 
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where maximum blade temperatures have been of the 
order of 1,300°F.‘ 


Precision-Cast Alloys 


The difficulty of forging the austenitic base alloys 
has led to the development of the precision-casting 
process for the preparation of turbine blades. Any 
large-scale production would so tax forging capacity 
that it would be almost mandatory to use cast blades. 
This casting process allows the use of many alloys with 
attractive properties for blades that are extremely dif- 
ficult or impossible to forge or whose good properties 
would be destroyed by forging. In many respects the 
properties of the cast form of alloys may be inherently 
more suitable for blade service than the wrought 
form. 

As in the case of forged alloys, it is difficult to gener- 
alize the metallurgical aspects of cast alloys. Such 
generalizations as are possible may be summarized as 
follows: 

(1) At temperatures somewhere between 1,350° 
and 1,500°F., the high-temperature strength of the 
cast form of most alloys becomes superior to the 
wrought. 

(2) The properties of cast materials may be in- 
fluenced by the pouring temperature, mold tempera- 
ture, and casting size and shape. 

(3) The metal and mold-casting temperatures 
control the grain size and carbide network in the cast- 
ing, as has been shown by Grant.> In some alloys 
the former controls the strength and the latter the 
ductility in the rupture test. 

(4) The grain size is also a function of the section 
cast, so that in a blade small grains exist in thin sec- 
tions and larger grains in the heavier section. 

(5) In large-sized blades difficulty may be encoun- 
tered in making sound castings and from large grain 
sizes. Large grains are believed to be associated 
with low fatigue strength, as well as other adverse 
mechanical properties. 

(6) When only a few grains occur in the casting, 
widely divergent properties may result, depending on 
the orientation of the grains. 

(7) Improvement in ‘strength may or may not 
result from solution treatments and aging. 

(8) Grant’s work has clearly shown the type and 
degree of control necessary for uniformly high prop- 
erties in one alloy. Such control is not easy in pro- 
duction without the added difficulty of grain-size varia- 
tion with section size. Limitations are also imposed 


by the need to fill both thin and heavy sections of a 
mold simultaneously and by the need for good surface. 
(9) The various alloys that have been developed 
from Vitallium show the effect of composition varia- 
tions much more clearly than do wrought materials. 
In experimental work with constant casting size and 
suitable control of the casting conditions, it has been 
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possible for Grant to show the effects of variations in 
chemical composition.® 

(10) It should be recognized that the cast alloys of 
the Vitallium type do not, in general, exhibit superior 
properties to the good wrought alloys until the tem- 
perature exceeds about 1,350°F. Their use at lower 
temperatures is, therefore, dictated by the advantages of 
the casting process and not by the properties of the 
alloys. 

Vitallium has been most widely used for actual 
blades. The reasons are that it was available at the 
time the need was greatest, and it has exhibited, in 
general, adequate properties. The other alloys de- 
scribed all have certain advantages and disadvantages, 
Some of them are now being used, and, as production 
difficulties are overcome and service requirements 
increase, it is certain that Vitallium will be replaced. 


FLAME TUBES 


The properties required of a metal for the flame 
tubes of combustion chambers are not clearly under- 
stood. This service is generally classed as an ex- 
tremely high-temperature low-stress application. 
Flame tubes currently have the shortest life of the 
high-temperature components of jet engines. The 
usual source of failure is direct flame impingement on 
the metal. It is extremely unlikely that any commer- 
cial metal would ever withstand the high temperatures 
and heat shock of direct contact with the flame. 

There are two principal reasons for flame impinge- 
ment: One is unsteady flame characteristics; the other 
is warpage of the metal into contact with the flame. 
The ultimate solution of the problem appears to be a 
matter of combustion engineering and design of flame 
tubes. 

The metallurgical characteristics involved are con- 
siderably different from those in the rotor. Insofar as 
it has been possible to analyze the service conditions, 
it appears that they consist of the following: 

(1) Adequate oxidation resistance for the normal 
operating conditions variously estimated from 1,650° 
to 1,900°F. 

(2) Good formability and weldability to permit fab- 
rication of the complicated shapes of flame tubes. 

(3) Resistance to warpage to prevent buckling and 
the formation of hot spots. This tendency is believed 
to be reduced by high thermal conductivity, emis- 
sivity, and perhaps high creep strength. 

(4) Other considerations may be resistance to heat 
shock without cracking and fatigue strength. 

The first two considerations have largely dictated 
the alloys used for flame tubes. These alloys are 
listed in Table 6. The high oxidation resistance and 
comparatively good fabrication characteristics of In- 
conel and 25Cr-20Ni + Si steels are well known. Vital- 
lium sheet on an experimental basis was successful, 
although the reason is not entirely clear. In England, 
Nimonic 75 has been used.” 
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TABLE 6 
Flame Tube Materials 
Chemical Composition, Per Cent ~ 
Alloy © Si Mn Cr Ni Co WwW Fe 
Inconel 0.15 0.50 1.0 13 75 > és 9 
25-20 + Si 0.05 2.0 1.8 25 20 he - 52 
Sheet Vitallium 0.20 0.6 1.0 25 3 61 5.5 2 
Nimonic 75 Not available 
TAIL PIPES TABLE 7 
5 In in 1,000-H Rupt i 
Metallurgically, materials for tail pipes and cones 
are similar to those for the combustion chamber, except 1,000-Hour Rupture 
that the temperatures are lower. The major problem Temperature Strength 
is the prevention of warpage during the rapid heating 198 + bl desis (°F.) (Lbs. per Sq.In.) 
and cooling during service. Oxidation resistance is Alloy and 19-9DL ,200 17,000 to 35,000-50,000 
14,000 to 20,000 


not so critical as in the flame tubes. For this reason, 
Type 347, Columbium stabilized stainless steel, has 
been widely used. 


Nozz_LE ASSEMBLY 


The requirements for nozzle guide vary to include 
oxidation, distortion, creep, and erosion resistance. 

The oxidation resistance of the materials used must 
be sufficient to prevent scaling at temperatures up to 
that of the gas stream. As low a coefficient of ex- 
pansion as possible and ability to withstand sudden 
temperature changes are needed for distortion resist- 
ance. Creep strength may be important in distortion 
resistance and certainly must be sufficient to withstand 
the gas bending stresses. Attention must also be 
given to the selection of an alloy that will not pit and 
erode under the action of the gas stream. Heat- 
shock resistance without cracking is also necessary. 

The service conditions are in some respects not so 
severe as for the rotor blades because the stresses are 
lower. The maintenance of throat area and finish of 
nozzle blades are important to engine performance, and 
for this reason consideration must be given to erosion 
and distortion resistance. 

In general, the same materials as are used for rotor 
blades would be satisfactory. Lower creep require- 
ments, however, have allowed the use of such alloys 
as 25Cr-20Ni alloy. Precision castings seem to be 
ideal for stator blades and are being used more widely. 


FUTURE DEVELOPMENTS 


Any prophecies about the materials that may be 
developed in the future are difficult to make. Some 
interesting developments, however, appear to be com- 
ing. These will not be purely metallurgical in nature. 
In fact, the major developments of the near future 
will probably be design improvements to take better 
advantage of present materials and to avoid their 
shortcomings. Such factors as combustion control, 
control over operating temperatures, and the use of 
cooling devices will serve to reduce the demands on 
the materials. 


Inconel W to Inconel X ,500 8,000 to 18,000 


1 
$495 to S590 and S816 1,500 

1 
Vitallium to X-40 1,500 14,000 to 25,000 


TABLE 8 . 
Selected Rupture Properties of Chromium-Iron-Molybdenum 
Alloys! 


Rupture Strengths 


Temperature (Lbs. per Sq.In.)——— 
Alloy °F.) 100 Hours 1,000 Hours 
60Cr-25Fe-15Mo 1,350 68,000 52,000 
1,600 20,000 
55Cr-25Fe-20Mo 1,350 73,000 54,500 
1,600 24,000 


One factor unique to the aircraft engine will prob- 
ably have important influence on the whole high- 
temperature metallurgical field. The essential re- 
quirement of saving weight will force better metal- 
lurgical control of high-temperature properties. At 
the present time, it is doubtful if there is sufficient 
knowledge of the variables involved or their effects 
on properties so that a producer can work to a narrow 
rupture-strength or creep-strength specification for the 
new high-strength alloys. This type of information 
should be obtained in order that production procedures 
can be developed which will provide the special prop- 
erties required for each application. The informa- 
tion is tedious and expensive to obtain and is even 
more difficult to apply in production. Yet it is a 
necessary requirement to efficient utilization of alloys. 
Efficient utilization, of course, means safe design with- 
out excessive safety factors and the use of the lowest 
cost materials possible. 

While the efficiency characteristics of the pure jet 
engine may not require higher and higher temperatures 
for better performance, the development of the gas 
turbine for driving propellers and other types of en- 
gines certainly will do so. Undoubtedly, there will be 
new alloy developments that will match those of the 
last few years in extending the permissible tempera- 
ture range—for example, those in Table 7. 

Entirely new base alloys will be developed to the 
point of useful application. The chromium-base al- 
loys, Table 8, pioneered by the Climax Molybdenum 
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Company, under the auspices of the O.S.R.D., is a 
good example of a new alloy system that should pro- 
vide materials for higher stresses or higher tempera- 
tures. The development of such materials is, however, 
a long and tedious process involving the gradual sur- 
mounting of melting, casting, and fabrication dif- 
ficulties in the laboratory and then in production. It 
seems to be characteristic of these new refractory alloys 
that much work must be done to develop toughness 
and that most things that help are detrimental to 
strength. 

Ceramic materials theoretically offer distinct ad- 
vantages for both pure ceramic construction and as 
coatings for metals. Progress to date has been en- 
couraging. Low heat-shock resistance and_ stress 
concentration brittleness are shortcomings that re- 
main to be conquered. There are also many problems 
in fabrication which will have to be solved. 


Perhaps the extension of service temperatures up- 
ward, beyond the limits of present alloy systems, will 
depend on the gradual improvement of the normal 
engineering characteristics of refractory metals and 
ceramics, combined with a new concept of design and 
fabrication to allow the use of brittle materials. 
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Preboiling Evaporation Losses in 


Aircraft Fuel Tanks 


L. M. WHITNEY, F. G. BOLLO, ano A. G. CATTANEO 
Shell Development Company 


SUMMARY 


In view of increasing interest in methods by which fuel evapora- 
tion losses can be predicted for long-range cruising flights at 
below-boiling altitudes, a set of equations was derived by which 
the losses at altitudes lower than the boiling altitude may be 
calculated from known initial fuel characteristics and known 
flight and tank conditions. 

The variables that affect these losses can be classified as fuel 
characteristics and temperature, altitude pressure, and 
tank and flight conditions. Some assumptions were necessary 
to derive the equations and to simplify them for practical 
use. 

There are three main causes of evaporation losses: tank 
breathing, air sweeping through the tank, and air evolving from 
the fuel. Equations were derived to determine the losses for 
each of these. Companion equations, needed to complete the 
solutions, were derived for the change in fuel vapor pressure re- 
sulting from fuel evaporation and its effect on temperature. The 
effect of evaporation on vapor pressure was determined experi- 
mentally, using a vapor pressure apparatus more sensitive than 
that used in the Reid method. For the temperature effect, an 
equation already determined by others was used. 

In practice, losses may vary in a single climb from less than 
1 to as high as 2 per cent, to which tank breathing contributes 
the greater part for tanks more than 15 per cent empty. If air 
sweeping is possible, the losses increase with the length of the 
flight. 


NOMENCLATURE 


P, = air-free partial pressure of gasoline, in. Hg. (Estimated 
from Reid vapor pressure and Fig. 2) 

P,; = total pressure (in. Hg) 

P, = partial pressure of air (in. Hg) 


K = air solubility coefficient (Ibs. of air dissolved per 100 
Ibs. of liquid fuel per in. Hg of air pressure) 

L_ = loss (per cent by weight) 

M, = molecular weight of gasoline vapor (assumed to be 
78.5) 


M, = molecular weight of air (28.9) 
W, = weight of liquid gasoline (Ibs.) 
Wa. = weight of air which becomes saturated with gasoline 


vapor (Ibs.) 

a = weight of air which becomes saturated with gasoline 
vapor per unit pressure. P will depend on the type 
of evaporation loss 

C = volume of fuel in tank (gal.) 

d = specific gravity 

F, = air flow (Ibs. per min.) 

T = temperature (° F. Abs.) 

h = time (min.) 

V = volume of vapor space in tank (gal.) 


Received February 7, 1947. 


INTRODUCTION 


EVAPORATION LOSSES from aircraft fuel 
tanks during flight have been given considerable 
attention. The work has been mostly concerned with 
losses occurring under conditions at which the fuel 
boils. Empirical equations have been obtained by 
which fuel losses are calculated from changes in fuel 
properties, such as specific gravity and tetraethyl lead 
content, or from flight altitude and fuel tempera- 
ture. 

Data obtained in studying evaporation losses after 
boiling show that some losses occur before boiling. As 
the following will show, these losses can be as high as 
2 per cent of the total fuel, depending on conditions. 
For a large aircraft carrying 3,000 gal. for a long-range 
flight below boiling altitude and temperature, this may 
mean a loss of 60 gal. or 360 Ibs. Considerable long- 
range cruising, especially over ocean areas, is done be- 
low the boiling altitude of the fuel, and a need exists 
for a method of calculating preboiling evaporation 
losses. 


DISCUSSION 


In general, these preboiling losses are due to the 
expulsion of fuel-saturated air from the tank vent, 
which may be the result of tank breathing, air sweep- 
ing, or air evolving from the fuel. As a step toward 
obtaining equations by which evaporation losses may 
be calculated, it is assumed that the losses are mainly 
due to these three causes. 

There are many variables that affect the evaporation 
loss of fuels in flight. These fall into four main cate- 
gories: (1) fuel characteristics, (2) fuel temperature, 
(3) altitude pressure, and (4) tank and flight condi- 
tions. In order to derive equations usable in practice, 
many of the variables whose effects are small need not 
be considered, and several assumptions can be made to 
simplify the equations. 


VARIABLES AFFECTING EVAPORATION LOSSES 


Fuel Characteristics 


Fuel characteristics that affect evaporation losses 
are vapor pressure and specific gravity of the fuel, 
saturation and solubility with air, and molecular weight 
of the gasoline vapor. 
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The molecular weight of the gasoline vapor was as- 
sumed to be 78.5. This figure is based on experimental 
work! that shows that the molecular weight of fuel 
vapor for current aviation gasolines remains constant 
up to 10 per cent evaporation losses. 

Vapor pressure of gasoline is usually measured by the 
Reid method. However, for work concerned with small 
variations in vapor pressure a method that will be more 
accurate than the present Reid vapor pressure method 
must be used.? Vapor pressure of fuels at different 
temperatures can be calculated by use of the slope of 
the A.S.T.M. distillation curve at the 10 per cent 
point.? 

For a given set of conditions losses will vary with the 
vapor pressure of the fuel. Fuels with low vapor pres- 
sures will have low losses. Evaporation of the fuel 
lowers the vapor pressure, which consequently helps 
to reduce the losses. 

Experimental work was necessary to obtain data for 
deriving a relationship between the change in vapor 
pressure and evaporation loss. An empirical equation 
was obtained for this relationship to be used with other 
equations for calculating losses. For small losses up 
to 0.4 per cent weight, the change in vapor pressure 
due to evaporation is sufficiently small so that it can 
be assumed that the vapor pressure remains constant. 
The resulting error was less than 0.1 per cent. 


Fuel Temperature 

A low initial fuel temperature decreases the preboil- 
ing evaporation losses and increases the altitude at 
which boiling losses occur. For example, a 6 lb. per 
sq.in. Reid vapor pressure fuel flown to 15,000 ft. 
loses 0.25 per cent by weight if its initial temperature 
is 90°F., compared with 0.45 percent at 100°F. and 
0.75 percent at 110°F. 

The initial fuel temperature is changed by heat loss 
to the walls of the tank and by the temperature drop 
due to evaporation losses. Actual data taken in flight 
with self-sealing tanks show drops in temperature of 
less than 1°F.* On the basis of calculations, heat 
losses to the walls of self-sealing tanks will drop tem- 
peratures on an average of 0.03° F. after 20 minutes’ 
flying time. The temperature drop ascribable to evapo- 
ration losses tends to decrease due to the heat retained 
by the walls of the fuel tank. This is especially true in 
rapid climb where the tank walls remain substantially 
at the initial temperature. 

Experimental work on drop in temperature caused 
by evaporation losses shows that the change in tem- 
perature, °F., is numerically equal to about 2.75 times 
the weight percentage lost.2 This temperature drop 
also reduces-the vapor pressure of the fuel. 


Altitude Pressure 

The loss of fuel increases with increasing altitude, 
since decreasing altitude pressure decreases the partial 
pressure of air in the vapor space of the fuel tank and 
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dissolved in the fuel. It is assumed that the pressure 
in the fuel tank is equal to atmospheric pressures, which 
is justifiable when considering rates of climb used in 
commercial operations. For example, laboratory work‘ 
shows that changes in the rate of climb up to 2,000 ft. 
per min. changed the evaporation loss less than 0.1 
per cent weight of the total fuel for data taken up to 
20,000 ft. altitude. 


Tank and Flight Conditions 


Experimental work done in simulated flight tests 
shows that fuel depth and fuel surface area have little 
effect on evaporation loss.‘ 

There then remain three main causes for evaporation 
losses from fuel tanks: (1) tank breathing, which is 
the loss in vapors from the fuel tank caused by changes 
in pressure; (2) air sweeping through the tank; and 
(3) air evolving from the liquid gasoline. It is as- 
sumed that the dissolved air comes out in equilibrium 
with the altitude pressure because of agitation by the 
booster pump. If, in the case of supersaturated fuel, 
the air does not evolve at the equilibrium altitude, 
there will be a greater loss as more air evolves at a 
lower altitude pressure. 

It is assumed that the rate of air sweeping the tank 
is constant regardless of altitude. For further simpli- 
fication it is assumed that the rate of climb is constant 
in terms of minutes per inch of mercury, which for alti- 
tudes below the boiling altitudes of the average aviation 
gasoline leads to very small errors. Losses from air 
evolving will occur in conjunction with both tank 
breathing and air sweeping the tank. 


GENERAL DERIVATION OF EQUATIONS 


After making the above assumptions, the evaporation 
loss can be expressed as a function of weight of air 
saturated with gasoline vapor per unit pressure, fuel 
air ratio, and change in partial pressure of air. 


L= F{(W,/P), (M,/M,), (P,/Pa); 


The subscripts, g and a, of the symbols stand for 
gasoline and air; for example, P, is the partial pres- 
sure of air. 

The loss in weight of fuel, whether by tank breathing 
air sweeping, or air evolving, resulting from a small in- 
crease Ai of the altitude i can be expressed as follows: 

W,M,P, 


, 


The total loss of fuel as Ai approaches 0, going from 
initial altitude pressure P,, to altitude pressure P,,, 


Pa ai) 


W,M,P 
Wa — Wut = p. Pa at ay) = 
[(WeM,P, 
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PREBOILING EVAPORATION LOSSES IN FUEL TANKS 705 


The quantities M,/M, and W,/P are assumed con- 
stant, and for a particular set of conditions sufficient 
data will be available to evaluate them. The equation 
for total loss of fuel is then: 


W,M, 
Total loss of fu P M, P, af, 
Wa M, Px 
ru og P., 


The limits are initial altitude pressure, P,, to final 
altitude pressure, P,,. Since AP, is an infinitesimal of 
higher order than AP,, P, is considered constant for 
very small changes in P;. Therefore the relationship 
P, — P, may be substituted for P,. 


= P, F; 
W.M, Pp 

Total loss of fuel = > uM, P, log. (P: — Py) P,, 
General Equation 

Total loss of fuel = 

(Pe (Pu — Pas) Pry (Pry Pos) 
Using M, = 78.5, Mz = 28.9. Converting to logi 
and per cent loss the equation becomes 
75.10 Wa 


[Poo logio (Pio — Po.) — Pon logio (Pin — Pon)] (1) 


For small losses, less than 0.4 per cent, it is assumed 
that P,, = P,, and the equation becomes simplified 


Simplified Equation 
75.10 W, — Pru | 
L Cd P Pa logio P,, ( ) 


For small losses greater than 0.4 per cent, Eq. (2) is 
solved simultaneously with the equation derived for 
change in vapor pressure with loss which was obtained 
from data on weathered fuel samples by use of the 
vapor pressure apparatus. This equation is not ac- 
curate beyond 5 per cent (Fig. 1). 


Equation for Change of Vapor Pressure with Loss 


Pon = 0.97Pp — 0.2L + 0.08 (3) 


There is also the change in P,, due to change in tem- 
perature with evaporation loss. The equation for 
change in temperature with loss follows. 


Equation for Change in Temperature with Loss 
AT = 2.75L + 0.014L? (4) 


The change in P,, due to AT is obtained by using 
Fig. 2, which gives the change in vapor pressure with 
temperature.” 
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EVAPORATION LOSS, L, 
PER CENT WEIGHT 


Fic.1. Change in vapor pressure of aviation gasolines with small 
evaporation losses. 


SPeciIFIC EQUATIONS 


The general Eq. (1) or simplified Eq. (2) can be used 
to calculate the three main types of preboiling fuel 
losses by use of the value of W,/P applicable to the 
particular situation and simplified by use of the as- 
sumptions discussed. Specific equations as obtained 
for each case are given below. 


A. Tank Breathing 


W,/P = initial weight of air in vapor space of 
tank (Ibs.) divided by partial pressure of 
air at sea level (in. Hg) 


397.47V Pu 
REID, PSI 8 
30 7 
© 
x 6 
Say 
$20 
a 
A 
| 
«x 10 4 
ro) MAW A 
< 
me 


60 80 100 120 140 160 
TEMPERATURE, °F 


Fic. 2. Change in vapor pressure of aviation gasolines with 
temperature. (Data calculated from reference 2.) 
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P,, = partial pressure of air at standard pressure 
= (29.92 — 


(29.02 — dTa = Pye 
(5) 


B. Air Sweeping 


As previously discussed, the rate of air sweeping the 


_ tank and the climb expressed as time per unit altitude 


pressure are considered constant. 

W./P = weight of air sweeping tank (Ibs. per min.) 
multiplied by rate of climb (time per unit altitude pres- 
sure, in. Hg) = — 


75.10 h Py - 
= F P 6 
For airflow at level flight, P, is constant so that the 


loss of fuel can be expressed as follows: 
Total loss = F,h(M,/M,)(P,/Pa) 


(6a) 
L = (32.6/Cd) [Pon/(Pim — Pon)] 
This shows that the loss at level flight increases rapidly 
with h.* The value of P,, is obtained from Eq. (3). 
For large values of h the method is not applicable due 
to the means by which Eq. (3) was derived. 


C. Air Evolving 


W,/P depends on solubility constant K for air dis- 
solved in gasoline given in terms of weight of air per 100 
weight of liquid fuel per inch of Hg air pressure.” 


W,/P = K = 0.00695 (1 — 1.1254) 


L = 0.043 (1 — 1.125d) | logio 
(7) 


To determine the weight percentage of fuel lost, the 
following data are needed: 

Fuel Characteristics—Reid vapor pressure, lbs. per 
sq.in.; specific gravity; initial fuel temperature, ° F.; 
slope of A.S.T.M. distillation curve at 10 per cent evapo- 
rated in ° F. per cent = Ti5y%, — T5q,/10; 

Tank and Flight Conditions.—Tank breathing, in- 
cluding capacity of tank, amount of fuel in tank, ini- 


tial and final altitude of flight; air sweeping, including 


* Actually the design of current fuel systems does not permit 
air sweeping through the tank. 
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EVAPORATION LOSS, PER CENT WEIGHT 


Fic. 3. Calculated preboiling evaporation losses of aviation 
gasoline. 


amount of fuel in tank, airflow, rate of climb, initial 
and final altitude of flight; air evolving, covering ini- 
tial and final altitude of flight. 

Preboiling evaporation losses can be calculated by 
placing the values obtainable from the above data in 
Eqs. (5) through (7) and solving simultaneously with 
Egs. (3) and (4). 


CORRELATIONS 


Experimental data on preboiling losses, obtained by 
others under simulated flight conditions incidental toa 
study of evaporation losses after the fuel boils, correlate 
well with values calculated by the above method. In 
one case, air evolving from fuel up to 14,500 ft. gave 
0.06 per cent losses compared to calculated values of 
0.07 per cent. For air breathing, losses of nearly 1 
per cent by weight were found on flights at 2,000 ft. 
per min. climb up to 16,500 ft.,* compared with calcu- 
lated values of 1.2 per cent. 

Calculated losses at different altitudes for fuels hay- 
ing 6, 6.5, and 7 Ibs. per sq.in. Reid vapor pressures are 
shown in Fig. 3. 
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Stress Analysis of Noncircular Rings for 
Monocoque Fuselages’ 


CHI-TEH WANG? ano S. RAMAMRITHAM? 


New York 


SUMMARY 


The application of the method of least work to the analysis 
of fuselage rings admits two possible sources of computational 
error: the errors introduced in solving the simultaneous equa- 
tions and the errors due to the final expressions for the bending 
moments around the ring being obtained, as a rule, in terms of 
differences of large numbers. These errors are generally recog- 
nized, and methods* ® 6 are available either to eliminate the first 
source of error or to reduce the second source. A method that 
reduces both sources of error is suggested in the present paper, 
and the advantage of such a procedure is illustrated by a numeri- 
cal example. 


INTRODUCTION 


a STRESS ANALYSIS OF RINGS for monocoque fuse- 
lages has a considerable amount of literature. A 
closed ring loaded in its own plane is statically inde- 
terminate to the third degree. The classical procedure 
in the analysis of such an indeterminate frame has been 
to make an imaginary cut through the structure, 
thereby reducing it to a statically determinate system. 
The three unknowns—viz., the bending moment, the 
axial force, and the shear force at the cut—can be cal- 
culated by strain energy methods. Such a procedure 
for analyzing rings of arbitrary shape is described by 
Lundquist and Burke! and has been discussed in detail 
by many other writers. The method consists of setting 
the partial derivative of the total strain energy of ring 
with respect to the unknowns equal to zero, according 
to Castigliano’s theorem, and then solving the resulting 
simultaneous equations for the unknowns. 

There are two possible sources of error involved in 
these computations. One is the errorinsolvingthesimul- 
taneous equations, where a high standard of accuracy 
is necessary to eliminate cumulative errors.? The 
other is due to the fact that the final expressions for the 
moments about the ring are given, as a rule, in terms of 
differences of large numbers; therefore, the unknowns 
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Unwersity 


cannot be obtained accurately by using merely the 
significant figures of a slide rule.* It was actually sug- 
gested by Ruffner‘ that all calculations should be car- 
ried out to five significant figures in the case of circular 
rings. But for noncircular rings, where it is usually not 
possible to obtain analytical expressions for the shape 
of the ring and the distribution of moments, the use of 
five significant figures becomes meaningless. 

In order to reduce the amount of work involved in 
the computations of the unknowns and also to reduce 
the first source of error, attempts have been made by 
Liepmann® and DuPlantier® to obtain the unknowns 
for a noncircular ring directly by corresponding indi- 
vidual expressions, thereby avoiding the solution of 
simultaneous equations. Liepmann follows the Miiller- 
Breslau’s procedure and assumes the unknown quanti- 
ties to act at the “‘elastic center’’ of the ring, and, from 
the condition that continuity must be reproduced 
at the point of cut, he obtains direct expressions for the 
unknowns. DuPlantier suggests the use of Hardy 
Cross’s column analogy method, where the structure is 
made determinate by a single cut, and the indetermi- 
nate moments are obtained directly by picturing an 
analogous column loaded eccentrically. The moments 
in the simple structure are obtained from the principles 
of statics. The difference of the two moments is the 
final required moment around the ring. The normal 
and shearing forces in the ring can also be easily ob- 
tained by simple calculation. 

To improve the accuracy due to the second source of 
error, Hoff* suggests reducing the ring to a statically 
determinate structure by introducing imaginary hinges 
at three points along the ring, instead of making a single 
complete cut. If the three imaginary hinges are located 
at the points where the final moments are small, then 
the final expression for the moment is the difference of 
one large value and one smallone. The improvement in 
accuracy by the above device isevident. However, in 
Hoff’s treatment, the three indeterminate moments at 
these imaginary hinges are calculated by the method of 
least work, as is done by Lundquist and Burke. In this 
way, not only is the amount of work increased, but the 
first source of error is also introduced. 

A method is suggested in this paper to reduce the ring 
to a statically determinate structure by inserting three 
imaginary hinges near the inflection points, as proposed 
by Hoff, and to determine the indeterminate moments 
by the column analogy method. It is evident that the 
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present approach possesses all the advantages of Hoff’s 
method and the column analogy, in that it reduces both 
sources of computational error, permits rapid semiauto- 
matic computation free from any sign confusion, and is 
applicable to rings of any arbitrary shape. 

It is assumed throughout this paper that the elastic 
limit of the material is not exceeded, that the strain 
energy due to axial and to shear forces is negligible in 
comparison with that due to bending moment, and that 


_ the direct stress distribution over the cross section due 


to bending is linear. Bending moments are considered 
positive if they cause compression on the inside of the 
ring. 


AX 
7 \ i 
7 
/, 
A' 
/ - 
AY, 
B 
Fic. 1. 
THEORY* 


Consider a closed ring with the Young’s Modulus £, 
the moment of inertia J, and the curvature either con- 
stant or varying as shown in Fig. 1. Suppose a bending 
moment is applied at the section B. If now the con- 
tinuity of the ring is destroyed by cutting it at the point 
A, the chord BA will swing through the small angle 
dd, the upper side of the section A moving to A’, while 
the lower part of the ring is assumed fixed. Since the 
angular movement caused by loads is small, the line 
AA’ may be considered as perpendicular to chord BA. 
From similar triangles, 


Rotation = do 
Vertical displacement A (do) 
1 e = — = 
sp V1 BA™ (1) 
Horizontal displacement Ax; = = (dd) 


But if we consider a more general loading condition, 


* A brief résumé of the method is hereby included for reader’s 
convenience. Detailed discussions can be found in DuPlantier’s 
paper® or any standard text on statically indeterminate struc- 
tures such as reference 2. The sign convention used in these 
references, however, is not correct, and this fact has been pointed 
out by Hoff.’ 


say, a system of coplanar forces acting on the ring, and 
if M, represents the moment at any point in the simple 
structure due to the external loading, we can obviously 
obtain the total rotation @ by using the equation 1/R = 
d0/ds = M,/EI as follows: 


Sf (M,/EDds (2) 


where R is the radius of curvature. Assuming a system 
of coordinate axes with origin at A with using Eq. (2), 


Total rotation = Sf (M,/EDds 
M,xds 
1 ical di = 
Total vertical displacement J EI 03) 


M, d. 
Total horizontal displacement = J 


EI 
Again, if M, represents the indeterminate moment at 
any point (x, y) on the ring due to the redundant un- 
known forces acting at the cut, then the total rota- 
tional and translational displacements at the cut will 
obviously be given by 


Total rotation = —Sf 
Total vertical displacement = — J ae : (a 


M, wds 
EI 


In order that continuity may be restored at the cut, the 
displacements represented by Eqs. (3) and (4) must be 
equal and, hence, 


JS (M.ds/EI) + = 0 
+ [Gry 
S Car) + {Gr 
Replacing the terms M,ds/EI and Mjds/EI by the 


symbols @, and 6,, respectively, and the integrals by 
summation signs, the Eqs. (5) can be rewritten as 


=O, +26; 


Total horizontal displacement = — J 


=0 (6) 


+ =0 


If we think of the @ values as vertical forces acting on 
any otherwise unsupported body, then the Egs. (6) will 
express the fact that this body is in equilibrium under 
the effect of these @ forces. 


The Analogy 


The statically determinate @, values may be looked 
upon as loads on the analogous column (Fig. 2), whose 
cross-sectional area is 2 ds/EI; this cross section has 
an axis of length fds as in the original ring, but its 
width must be 1/EIJ at each point. The height of the 
column is immaterial. The values of 6; represent reac- 
tions for the analogous column. The tilting moments 
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caused by eccentricity of @, loads are = 6,x about the 
yy axis and 2 6,y about the xx axis. Finally, the values 
of M, = 6,/(ds/EI) represent unit stresses in the analo- 
gous column. The analogy is now complete, and the 
obvious conclusion can be drawn that the indeterminate 
moments M; can be calculated from the simple equa- 
tion for an eccentrically loaded column: 


W May . Myx 
where 
f =M; W=280,; A ds/EI 
My, = 2 6y; My = 6x 
ds ds 
Tes zy (=) 


where x and y are coordinates along principal axes 
through the centroid of the cross section of the analo- 
gous column. 


Sign Convention 


When the indeterminate moments are known, the 
final moments can be obtained by combining them with 
the statically determined moments in the cut ring— 
i.e., 


M=M,-™, (8) 


Signs follow from an easily remembered rule. Consider 
compression stresses in the column as positive, and add 
or subtract the terms of the column analogy Eq. (7) 
just as in an eccentrically loaded column—i.e., add 
compression stresses together but subtract any tension 
stresses. The total moment at any point is equal to 
the simple structure moment M, minus the indeter- 
minate moment 


PROPOSED METHOD OF ANALYSIS 


The proposed method for the stress analysis of non- 
circular fuselage rings for monocoque fuselages consists 
of the following steps. 

(1) Choice of the Simple Structure.—Referring to Eq. 
(8), it is shown that the final moments are the difference 
of M,and M,. The designer can visualize fairly clearly 
the deformations and deflections that would be caused 
by the external loads acting on the ring. From the de- 
flection curves a rough mental picture will give the ap- 
proximate locations of the points of inflection where the 
final moments will be small. By inserting the three 
imaginary hinges near the inflection points, the mo- 
ments M, in the simple structure are made close to the 
final values. Thus, the indeterminate moments M;, to 


be determined from column analogy are small compared 
to M,. 

(2) Determinate Moment M,.—Each of the three sec- 
tions into which the ring has been divided can be con- 
sidered independently as a system in equilibrium, and 
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the reactions at the pin joints can be computed from 
the three equations of equilibrium. In the case of 
monocoque fuselage rings, where the reactions are pro- 
vided by the tangential components of the shear forces 
in the skin, it is sometimes simpler to use the graphical 
method of finding the M, curve, as shown in reference 3. 

(3) Indeterminate Moments M;,.—This part of the 
final moment due to the continuity of the structure is 
obtained by the application of the column analogy 
method, and it consists of the following steps: 

(a) Geometric properties of the analogous column: 
The cross section of the analogous column is of such a 
shape that its centerline is identical to that of the ring 
under investigation, and its transverse width at any 
point is numerically equal to the reciprocal of the 
product, EJ, of the redundant structure at the corre- 
sponding point. Hence, the cross-sectional area of 
the analogous column is 


A=fdA=/YS (ds)(1/ED) 


Then the position of the centroid of the analogous 
column is found. By choosing a convenient pair of 
rectangular coordinate axes with the centroid as origin, 
the moments of inertia and the product of inertia of 
the column section are given by 


= J = JS (ds/Ely’ 
In = = S (ds/EDx? (9) 
= S xydA = JS (ds/El)xy 


(b) Total 6, loads on the analogous column: With 
the M, and the area of the analogous section known at 
every point, the total @, loads can be computed by 
analytical or graphical methods. 

(c) Tilting moments about the centroidal axes: 
The tilting moments M,, and M,, of the @, loads about 
the centroidal axes x and y, respectively, are given by 


Miz = S 6.9 = S(M.ds/El)y 


My, = S 0x = S (Mds/EDx (10) 
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FIGURE 5. 


(d) The indeterminate moment M,: Then the M, 
moments are given by the column analogy equation, 
since 


wi 
( gy ( ) 


f= W/A + Cu + Cy 


In the case of rings with an axis of symmetry, the prod- 
uct of inertia 7,, = 0, about the principal axes, and 
Eq. (11) is reduced to 
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The correspondence in the analogy has been explained 
previously, and the fiber stress in the analogous column 
is the same as the W, value in the ring. 

(4) Final Moment M in the Ring.—Thefinal moments 
are obtained by subtracting the indeterminate moment 
M, from the simple structure moment M,. It must be 
noted that the sign of the indeterminate moment is 
automatically determined by using the sign conventions 
that moments are positive if they cause compression 
on the inside of the ring and that the compressive loads 
on the analogous columns are positive. Depending on 
whether M;, is of the same sign or of opposite sign to 
that of M,, the final moments will be less or greater 
than the simple structure moment. 


NUMERICAL EXAMPLE 


In order to illustrate the application of the proposed 
method for analyzing rings of arbitrary shapes, the 
same problem worked out in DuPlantier’s paper is 
chosen. This will enable a direct comparison of the 
results by the two methods and bring out clearly the 
advantages of the present method. The problem is to 
compute the bending moment along the trapezoidal 
ring shaped and loaded as shown in Fig. 3 and having a 
constant cross section and a constant modulus of elas- 
ticity. 

Reactions in the simple structure: The trapeziodal 
ring is made statically determinate by inserting imagi- 
nary hinges at EZ, F, and G (see Fig. 4). Let Vi, Vo, Vs 
and H,, He, H3 be the vertical and horizontal reactions, 
respectively, at the hinges E, F, and G indicated in the 
diagram. 

From Fig. 5a, equilibrium equations give 


= V2; A, = A; 20H, = 2V; (13) 


From Fig. 5b, equilibrium equations give 


Fn = V2 Vs = 15H 20V2 = —10 (14) 
Substituting for H2, V2, and V; from Eq. (13), 
Vs; = 85/185 Ib. 


Section properties of the analogous column: Divide 
the ring axis into eight segments 10 in. long (As = 10 
in.) and two segments 14.14 in. long (As = 14.14 in.). 
Locate the centers of these segments and assign to them 
station numbers from 1 to 10 (see Fig. 3). These sta- 
tions apply also to the section of the analogous column. 
For each of these divisions compute their respective 
AA values, where AA = As/EI, computing the J value 
at each point and assuming it to be constant throughout 
the As length. In this problem both £ and J are con- 
stant and are assigned relative values of unity; hence, 
in column 2 of Table 1, AA = As. 

Choose any convenient pair of rectangular reference 
axes, x’-x’ and y’-y’ asin Fig.3. Place in columns 3 and 
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1 2 3 4 5 6 7 8 9 10 11 12 13 
Ring AA = z= Alz = zz: = Alyy AW x 185 
Stations As/EI x’ y’ AA (x’) AA(y’) AA (ry) AA (y?) AA (x?) M,X 185 (12) X (2) 
1,R 10 15 20 150 200 —2.08 11.85 — 246 1,404 43 — 800 — 8,000 
2 10 5 20 50 200 —12.0 11.85 —1,431 1,404 1,459 +200 +2,000 
3 10 0 15 0 150 —17.08 6.85 —1,170 469 2,917 +650 +6,500 
4 10 0 5 0 50 —17.08 —3.15 538 2,917 +550 5,500 
5, F 10 5 0 50 0 —12.08 —8.15 985 664 1,459 0 0 
6, S 16 15 0 150 0 —2.08 —8.15 170 664 43 — 1,000 — 10,000 
7 10 25 0 250 0 7.92 —8.15 — 645 664 627 —150 ,500 
8 10 35 0 350 0 17.92 —-8.15 — 1,460 664 3,211 +700 +7,000 
9 14.14 35 5 495 71 17.92 —3.15 —798 140 4,541 +750 10,600 
10,G 14.14 25 15 354 212 7.92 6.85 767 663 887 0 
TABLE 1 (Continued) 
14 . 15 16 17 18 19 20 21 22 23 24 
Moments by 
Ww -—Moments by Present Method—~ ——DuPlantier’s Analysis— 
AMyy X 185 AMzz X 185 za Mi = (16) + = M = 
(7) X (13) (8) X (13) = AA Cz(z) Cy(y) (17) + (18) M, M,- Mi M; M, M,—-Mi 
0.604 +0.088 —0.208 +0.484 —2.03 —2.51 2.54 0 —2.54 
16.64 X 103 —94.8 X 108 0.604 —0.063 —0.208 +0.333 —4.33 —4.66 4.68 0 —4.68 
0 0 0.604 —0.302 —0.208 +0.094 —0.094 
—24.16 X 10% 23.7 X 108 0.604 —0.357 —0.208 +0.039 +1.08 +1.04 8.98 10 1.02 
0.604 —0.516 —0.208 —0.120 3.78 +3.90 11.12 15 3.88 
—111 X 103 44.5 X 10° 0.604 —0.516 —0.120 —0.032 3.52 +3.55 11.47 15 3.53 
—94 X 10? —17.3 X 108 0.604 —0.516 +0.055 +0.143 2.97 2.83 12.18 15 2.82 
0.604 —0.516 +0.143 +0.231 2.70 +2.47 12.54 15 2.46 
0 0 0.604 —0.357 +0.143 +0.390 —0.39 10.40 10 —0.40 
20.8 X 103 81.5 X 103 0.604 —0.063 +0.143 +0.684 —5.40 —6.09 6.10 0 —6.10 
—11.88 X 10% 12.2 X 1 0.604 +0.239 +0.143 +0.986 —0.81 —1.80 1.81 0 —1.81 
125.4 X 193 —57.05 X 10% 0.604 +0.542 +0.143 +1.29 3.79 +2.50 —2.48 0 2.48 
0.604 +0.693 +0.143 +1.44 6.08 +4.64 —4.63 0 4.63 
190 X 103 —33.41 X 10% 0.6 +0.542 +0.055 +1.20 4.05 +2.85 —2.84 0 2.84 
0 0.604 +0.239 —0.120 +0.723 0 —0.72 +0.74 0 —0.74 


4 the x’ and y’ coordinates of the station centers, calling 
these values positive when to the right and above the 
reference axes. The centroid of the analogous column 
is located by applying the equations, 


2(AA)x’ column (5) 1,849 
= 17.08 in. 
SAA column (2) 108.3 
_2(AA)y’ column (6) 8.15 in 


and by extending columns 5 and 6. The centroidal 
axes x-x and y-y parallel to the existing axes are drawn, 
and the coordinates (x, y) of the stations are entered in 
columns 7 and 8, of Table 1. Using the values in Table 
1, the following properties of the analogous column 
section are calculated: 


Area of section A = 2AA =2(2) = 108.28 
Product of inertia J,, = 2(AA)xy = 2(9) = —3,290 
Moment of inertia J,, = 2(AA)y? = 2(10) = 6,835 


Moment of inertia J,, = 2(AA)x? = 2(11) = 18,104 


Moments in the simple structure (M,): With the 
reactions at the hinges known, the applied moment M, 
in the base structure may be computed, and entered 
in column 12. A positive value is assigned to any 
moment causing compression on the inside of the 
ring. 

Total 6, loads on the analogous column: The values 
of the product of VM, and the AA of the section stations 


along the ring giving the elemental 0, load on the section 
are entered in column 13. The sum of column 13 gives 
the total load on the column as 


W = M,(AA) = 12,100/185 
and 
W/A = (12,100/185)(1/108:3) = 0.604. 


The coefficients C, and C, of the column analogy 
equation: Referring to Eq. (11), 


C, (Teel yy — I,,*) 


where the tilting moment 
1,000 
(—94.8 + 23.7 + 44.5 — 17.3 4+ 


185 
81.5 + 12.23 — 57.05 — 33.41) 
— 220 F 


Il 


Mizz = 
The tilting moment 


My = =(M,)(SA)x 


(16.64 24.16 — 111 — 94.0 + 20.8 — 
11.88 + 125.4 + 190) 
M,, = 604 


By simple substitution, we have, 
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2 € 1 
3 10 / 
12°5 
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Ls 6 7 8 
D Ms (DUPLANTIER) 
75 
nK Ms (PRESENT METHOD) 
50 > L 
/ 
4 \ 
7 
\ Ve FINAL MOMENTS 
\} 
-5 VA 
RING STATIONS 
Fic. 6. 
604 X 6,835 — 220 X 3,290 — But, like any other mechanical device, such as the slide 
7 “~~ (6,840)(18,100) — (3,290)? or rule, there is considerable effort required for its mastery. 
Also, like the operation of the slide rule, the column 
pee (—220)(18,100) — (604)(—3,290) = 0.0175 analogy procedure, once learned, is easily retained in 
mind. At the same time it must also be remembered 


(6,840) (18,100) — (3,290)? 
The indeterminate moments M,: 
M, = W/A+ Cx + Cy 


The values of the first, second, and third terms of the 
right-hand side of this equation for the stations on the 
ring are entered in columns 16, 17, and 18 of Table 1. 
The values of the M, at every station are entered in 
column 19. 

Final moments: The indeterminate moment M/, in 
column 19 is subtracted from M, in column 20 to give 
the final moments, which are tabulated in column 21, 
and the results in Table 1 are plotted in Fig. 6. 


CONCLUDING REMARKS 


The column analogy method used in the present paper 
is a highly valuable tool for analysis of single-span 
structures with not more than three redundant elements. 
It should be looked upon as a sort of mechanical device 
that will provide a scheme with ease and directness and 
shorten the time needed for making computations. 


that all mechanical tools are dangerous in that they 
tend to relieve the computer of the need for maintain- 
ing what might almost be termed physical contact 
with the structure itself. 
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Cardboard-Box Wing Structures 


F. R. SHANLEY* 


Consulting Structural Engineer 


INTRODUCTION 


| para OF THE MILITARY RESTRICTIONS which ex- 
isted during the war, some of the interesting struc- 
tural developments have not been reported. One of 
these is the use of a “cardboard-box”’ type of metal 
structure, consisting of two sheets joined by another 
corrugated sheet, as shown in Fig. 1. 

The first use of this idea on a Lockheed airplane was 
made in 1940 on the P-38. The original wing structure 
of this airplane, in the fuel-tank area, consisted of a 
smooth outer sheet riveted to spanwise corrugations, 
with conventional chordwise ribs at fairly close inter- 
vals. When it became obvious that selfsealing fuel 
tanks must be installed for actual combat duty, the 
elimination of the ribs in the fuel tank area became 
necessary. 

The possibility of utilizing the chordwise stiffness of a 
panel to provide resistance to. buckling had been recog- 
nized,f but little experimental work had been done at 
that time. Rough calculations indicated that an inner 
sheet, riveted to the corrugations, might permit the 
ribs to be removed entirely from the tank area. The 
addition of the inner sheet completed the truss structure 
of a chordwise slice through the panel, thereby provid- 
ing chordwise bending rigidity. The additional sheet 
also increased the bending rigidity in the spanwise direc- 
tion and increased the amount of material in compres- 
sion. This, of course, reduced the working stress and 
made it possible to carry the required loads at a lower 
average allowable stress. 


TEST PANEL 


A test panel (Fig. 2) was made, simulating the com- 
plete panel over the fuel tank area. This was tested in 
end compression in a conventional testing machine 
with the edges of the panel supported as they would be 
in the airplane. The test showed that the required 
axial load could be carried by the panel without buck- 
ling. The final failure was local, not general. (See 
Fig. 2.) The ribs were therefore removed, and inner 
sheets were added on both top and bottom surfaces, 
using “‘blind’’ rivets. All P-38 airplanes were subse- 
quently made in this manner. (The material was 24S- 
T sheet and 24SR-T corrugations.) 


Received March 24, 1947. 

* Formerly Division Engineer in Charge of Engineering Re- 
search, Lockheed Aircraft Corporation. 

t See reference 1 for a discussion of the effects of chordwise 
rigidity. 
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EXPERIMENTS ON THE XP-58 


Because of the success of this experiment it was de- 
cided to use a similar structure on a new experimental 
airplane, the XP-58. This airplane, which did not go 
into production, was similar to the P-38 but consider- 
ably larger. (See Fig. 3.) To simplify production the 
wing was designed with a center section having constant 
chord, constant thickness, and no dihedral. The span 
of the center section was 36ft.8in. The primary struc- 
ture took the form of a box with three vertical spanwise 
shear beams. No internal ribs were used. The long, 


Fic. 1. Cross-section of typical sandwich panel. 


Fic. 2. Test panel for Lockheed P-38, after failure. Note local 
buckling. 
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Fic. 3. The Lockheed XP-58, which used the ‘‘cardboard-box”’ 
structure for the wing center section. 


Fic. 4. One cell of the box beam used on Model XP-58. 


narrow, selfsealing tanks were inserted from the twoends 
of the center section. Fig. 4 shows the general construc- 
tion of one cell, representing half of the box structure. 
The shear beams were made from flat sheet, stiffened by 
vertically corrugated sheet. 

The additional center shear beam was used because it 
was desirable to reduce the unsupported chordwise 
width of the panels. The corrugations otherwise 
would have had to be much deeper, would have become 
inefficient, and would have used up too much fuel-tank 
space. This, incidentally, is one of the design prin- 
ciples in the use of laterally supported panels. The 
important dimension is the width of the panel, which 
plays the same role that the rib spacing does for conven- 
tional construction. If the panel is kept relatively 
narrow and has chordwise bending rigidity roughly 
equal to, or greater than, its spanwise rigidity (per unit 
of length), there is nothing to be gained by reducing the 
unsupported length in the spanwise direction. This 
means that ribs are not required at all, so far as the com- 
pressive strength of the panel is concerned. 

On the XP-58 the nacelle and fuselage loads were 
introduced into the wing box structure by means of 
“external” ribs, which were simply a reinforced part of 
the nacelle and fuselage structure. These external 
structures also served to maintain the cross-sectional 
shape of the wing. In this hollow type of construction 
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it is necessary to have a minimum of two ribs or bulk- 
heads, one at or neareachend. They do not have to be 
internal, however. Where high local transverse loads 
are introtluced into the structure it is, of course, neces- 
sary to have ribs or their equivalent. 

In connection with the XP-58 design, a number of 
tests were made in the Lockheed Engineering Research 
Laboratory. These were for the purpose of determin- 
ing local buckling loads, transverse bending rigidity, 
general instability loads, and the effects of combined 
loadings (axial, shear, and normal pressure). Figs. 5 
and 6 show some typical test specimens. The theory of 
laterally supported panels was also worked out and 
checked by tests. Finally, a full-sized section of the 
XP-58 wing was tested under all critical loading con- 
ditions to ultimate loads. It is significant that the 
final failure of this test specimen occurred by local buck- 
ling near one end of the panel rather than by general 
instability. This showed that the absence of ribs was 
more than compensated for by the chordwise bending 
rigidity of the panels. 


DISCUSSION 


Although it is difficult to draw any general conclu- 
sions as to the strength/weight efficiency of the card- 
board-box type of construction, it is safe to say that, 
when properly used, it is at least no heavier than con- 
ventional construction. The average working stresses 
will probably be lower (an advantage from the stand- 
point of rigidity and fatigue). This does not mean that 


the weight is higher; it simply means that more material 
is put to work. 


In fact, the type of construction used 


Fic. 5. Test for transverse bending stiffness under normal 
pressure. (Lockheed Engineering Research Laboratory.) 


Fic. 6. A specimen that was tested to determine local buckling 
strength. (Lockheed Engineering Research Laboratory.) 
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in the XP-58 approaches 100 per cent utilization of the 
box-beam material in resisting bending and torsion. 

Although the greatest advantage comes in connection 
with the use of selfsealing fuel tanks, there are definite 
production advantages. On the XP-58 the top and bot- 
tom panels and the three vertical shear beams were 
fabricated separately, then joined together. The inner 
sheet was riveted to the corrugation by ordinary rivets, 
using a “magnetic mouse”’ type of bucking tool. (The 
tool was inserted into a loop of the corrugation, and the 
bucking force was obtained by a pull from a magnet at- 
tached to the rivet gun.) Incidentally, bullet pene- 
tration tests turned out quite satisfactorily for this type 
of structure. 


OTHER TYPES OF SANDWICH CONSTRUCTION 


In studies of various types of sandwich construction, 
several other ideas of a similar nature were tried out. 
One of these employed a Pratt truss arrangement of 
spacers and spanwise strips, as shown in Fig. 7. It 
was hoped that this arrangement would reduce vertical 
shear deflections, which had been found to be relatively 
important in the theory for general instability. Com- 
parative tests, however, showed that the corrugation 
type of sandwich would be more efficient. 

Another variation is shown in Fig. 8. Here the cor- 
rugation is replaced by spanwise spacers made from 
narrow bars. This type of construction does not have 
much chordwise rigidity with respect to transverse 
shear loads. It must therefore be used in narrower 
widths, or must have additional support from chord- 
wise ribs. It is also impractical to make unless the two 
skins are relatively close together. Nevertheless, 
there is a place for such sandwiches in areas where the 
unit axial loading (pounds per inch of width) is high. 
Incidentally, panels of this type, tested as simple col- 
umns such as shown in Fig. 8, showed that rivet spacing 
must be relatively close in order to make the two sheets 
fully effective as flanges of a beam. The spacing of the 
separators also affects the flange efficiency to some éx- 
tent; relatively wide spacing causes a considerable re- 
duction in efficiency. 

Many other types of metal sandwiches have been in- 
vestigated at Lockheed, including one in which two 
corrugated sheets were arranged to cross each other 
diagonally. For fairly heavy loading, however, none 
seemed to offer any advantages over the original type of 
cardboard-box construction. 


Fic. 7. Experimental panel (turned out to be less efficient than 
corrugated type). 


Fic. 8. A test panel employing narrow bars for spacers. (Lock- 
heed Engineering Research Laboratory.) 


CONCLUSIONS 


Perhaps the most important advantage of the true 
sandwich type of panel is its ability to resist bending in 
any transverse plane. It is more efficient to support 
such a panel along the edges than to break it up into 
short spanwise lengths by means of ribs. This requires 
a complete revison of basic design ideas, with ribs re- 
placed by spanwise “‘shear’’ beams or spacers. The 
same reasoning applies to the design of thick-skinned 
wings, in which the plate is heavy enough to be used 
without auxiliary stiffeners. In fact, the simple thick 
plate, supported by spanwise members of full wing 
depth, is the ideal type of wing construction from every 
point of view, except that for conventional airplanes the 
unit axial loading is too low to permit the use of thick 
plate at high allowable stresses. 

The idea of high allowable stresses must be used with 
caution, however, since it is not a direct measure of 
strength/weight efficiency. The only sure basis of 
comparison is a weight estimate of the entire wing struc- 
ture on a constant strength basis. (It may be incorrect 
to compare strengths on a constant weight basis.*) This 
might show that the wing having the lower allowable 
stress is actually lighter because more material is work- 
ing than in the conventional design. 

Finally, it should be mentioned that relatively high 
curvature (such as that of a leading edge) has roughly 
the same general effect as narrow panel width. Sand- 
wich construction can often be used for such elements 
without providing any further stiffening members (such 
as ribs). The development of lighter types of core 
material, now well under way, should ultimately permit 
the efficient use of sandwich construction in the fuse- 
lage, where the larger radius of curvature requires a 
large effective shell thickness in order to prevent buck- 
ling. 
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correspondence, avoid the return of papers for changes, minimize the work of preparation for the 
printer, and save the expense due to the charges made for ‘‘author’s corrections.” 


Manuscripts: The original typewritten copy of the paper is 
® desired, double or triple spaced on one side of white paper sheets, 

consecutively numbered. There should be wide margins to 
allow for the marking of directions to the printer. Correcting, 
changing, or adding to papers after they are intypeiscostly. It 
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Typographical errors may be corrected on proofs, but if authors 
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protection) to the Editorial Office, Institute of the Aeronautical 
Sciences, 2 East 64th Street, New York 21, N.Y. All manu- 
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initials of the author should be written as he prefers. The use 
of the full name of an author is advocated because of the fre- 
quent duplication of initials and surnames which sometimes 
makes it difficult to establish the identity of the author. This is 
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cases constitute the only source of information used in compiling 
scientific reference indexes. Abstracts printed in other journals, 
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B length of the paper. 
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tice. Too detailed tabular matter (general results of such tables 
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REFERENCES AND FooTNOTES: References should appear as 
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cannot be reproduced (including pencil drawings) will be returned 
to the author for redrawing, thus delaying publication of the 
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white contrasts. They must be on glossy white paper. Avoid 
round and oval photographs. 


CaPpTIONS AND LEGENDS: Legends or captions must accom- 
pany each drawing or photograph submitted. If written on the 
drawing or photograph, they should be placed below and well out- 
side the part to be reproduced. Each table should have a caption 
such as Table 1, Table 2, Table 3, etc. Captions should be com- 
plete in themselves so as to make the data intelligible to the 
reader without reference to the text. A duplicate list of captions 
for figures should be included as the last page of the manuscript. 
Use ‘Fig. 1” (not Figure 1), Figs. 3 and 4, etc., in both the text 
and the numbering of illustrations. In the text, ‘Eq. (1),” or 
‘*Eqs. (1) and (2)” are preferable to ‘‘Equation (1).” In captions 
and legends, except for ‘‘Fig.”” and “‘Eq.,”’ and in table headings, 
write all words in full; do not abbreviate. 


MATHEMATICAL WorRK: Only the simplest formulas should 
be typewritten; all others should be carefully written in pen 
and ink, the writing to be large enough so that ample room is 
provided to mark mathematical matter for the printer. A con- 
siderable space for marking should be allowed above and below 
all equations. All complicated equations should be repeated on 
separate sheets with plenty of space left for marking. The solidus 
should be used for simple fractions appearing within the text. 
Make all expressions clear to the typesetter. Greek letters used 
in formulas should be clearly designated by name on the margin 
of the manuscript. All symbols should be clearly written and 
carefully checked. The difference between capital and lower- 
case letters should be clearly distinguished and care taken to avoid 
confusion between zero (0) and the letter (0), between the numeral 
(one) and the letter (ell) and the prime (’), between alpha and a, 
kappa and k, u and mu, v and nu, nand eta. All subscripts and 
exponents should be clearly marked and dots and bars over letters 
or mathematical expressions should be avoided. Avoid compli- 
cated exponents and subscripts. When it is necessary to repeat a 
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tion in the 
nology which are applied to aeronautics. 


: Publications 


Papers presented at meetings or submitted to the Edi- 
torial Board are published in the JournaL or THE AERO- 
NauTicai Sciences or the AkRONAUTICAL ENGINEERING 
Review. The Journat prints in each issue several full- 
length scientific papers on new research and developments 
in various fjelds applied to aviation. The Review keeps 
members and subscribers up to date on aeronautical news 
and literature through reviews of new books and periodical 
articles, government publications and trade literature and 
— papers on applied engineering and aircraft pro- 

uction. 

The AgronavTicaL Enoingerinc Catatoc, published 
yearly, is a guide to sources and specifications of materials, 
parts, and accessories used in the design and production of 
aircraft, aircraft engines, and parts. 


Research Facilities 


Through the Aeronautical Archives of the Institute, 
members have for their use the most complete aeronautical 
library and reference research facilities maintained by an 
technical society. At the New York building of the Insti- 
tute, The W. A. M. Burden Library of books, reports and 
periodicals, the Aeronautical Index of subject files, bio- 
graphic files and bibliography, and the historical and art 
collections in the Archives are available for study. The 
Pacific Aeronautical Library at 6715 Hollywood Boulevard, 
Los Angeles, maintains a reference collection and reading 
room for members and loans books to company libraries 
in the area. Through The Paul Kollsman Lending Li- 
brary, aeronautical books are lent by mail without charge 
to members anywhere in the continental United States. 
In addition to a complete collection of standard aero- 
nautical reference works, the Kollsman Lending Library 
has available for loan new aviation books as they are 
reviewed in the AERONAUTICAL ENGINEERING Review. 


Membership 


In addition to membership in professional societies 
representing their roman fields, specialists applying 
their knowledge and experience to aeronautics find mem- 
bership in the Institute of great value in the contacts made 
with others engaged in aeronautical work and in the broad 
coverage of all aeronautical problems furnished through its 
meetings and ees: All applications for member- 
ship are carefully reviewed by the Membership Committee, 
which recommends the grade of membership for whicn the 
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nited States which brings together all of the sciences and branches of engineering and tech- 


applicant is found eligible. There are advanced grades of 
membership—Fellow and Honorary Fellow—to which 
members may be elected as their experience warrants. 
Thus, a member receives due recognition—more valuable 
because it is given by his colleagues—for his contributions 
to the progress of aeronautics. New applicants may be 
elected by the Membership Committee to the grade of 
Associate Fellow if they have been engaged for not less 
than ten years in the practice of the aeronautical sciences 
and shall have been in responsible charge of important 
scientific or engineering work, or shall have otherwise 
made outstanding contributions to the aeronautical 
sciences; to the grade of MEMBER if they have been en- 
gaged in professional work for at least eight years and have 
acquired a recognized standing in engineering, design, or 
other special work a to aeronautics; to the grade of 
Associate Member if they have acquired a recognized 
standing in an administrative capacity in the aviation in- 
dustry; to the grade of Technical Member if they are recent 
engineering school graduates or are engaged in technical 
aeronautical work. Student Members are admitted by 
application through Student Branches organized at their 
schools. Aeronautical companies, as well as individuals, 
are affiliated with, and participate in, the support of the 
Institute through Corporate Membership. 


Dues 


An entrance fee of $10 is required of all new members 
who apply for admission to a grade of membership to be 
specified by the Membership Committee (Corporate Mem- 
bers, at Members, and Student Members being ex- 
cepted ). 

Giant membership dues include a subscription to the 
AgronavTicaL ENGINEERING Review. Dues for the 
various grades of membership are as follows: Fellows— 
$16.50; Associate Fellows—$12.50; MEMBERS—$10.50; 
Associate Members—$10.50; Technical Members—(over 
26 years of age) $10.50, (under 26) $6.50. 

e regular subscription price for the JourNaL OF THE 
AERONAUTICAL Sciences and the AERONAUTICAL ENGINEER- 
1nG Review is $7.00 and $3.00 per year, respectively. 
Members may subscribe to the Journat at one-half the 
regular subscription rate. 

Those who are engaged in aeronautical work or who 
have an interest in any technical phase of aviation will find 
the services of the Institute to be a necessary aid in keeping 
abreast of developments in these times of such rapid prog- 
ress. An application form and further information about 
membership can be obtained from an Institute member or 
by writing to the Secretary. 
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